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Abstract In this paper, uniform pointwise regularity estimates for the solutions of conductivity 
equations are obtained in a unit conductivity medium reinforced by a e-periodic lattice of highly 
conducting thin rods. The estimates are derived only at a distance e 1+r (for some r > 0) away from 
the fibres. This distance constraint is rather sharp since the gradients of the solutions are shown 
to be unbounded locally in L p as soon as p > 2. One key ingredient is the derivation in dimension 
two of regularity estimates to the solutions of the equations deduced from a Fourier series expansion 
with respect to the fibres direction, and weighted by the high-contrast conductivity. The dependence 
on powers of e of these two-dimensional estimates is shown to be sharp. The initial motivation for 
this work comes from imaging, and enhanced resolution phenomena observed experimentally in the 
presence of micro-structures [53] . We use these regularity estimates to characterize the signature of 
low volume fraction heterogeneities in the fibred reinforced medium assuming that the heterogeneities 
stay at a distance e 1+T away from the fibres. 

Keywords homogenization - high conductivity - fibred media - weighted second-order elliptic 
equations - regularity estimates 

AMS subject classification: 35J15 - 35B27 - 35B65 



1 Introduction 

Consider a material contained in Q = u) x (— L, L) C K 3 where w is a bounded domain in K 2 with 
smooth boundary dm. Given some fixed loq <g ui, we assume that inside J?o = w o x {~L,L), the 
material contains small cylindrical rods of high conductivity. A sketch of the domain is represented 
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in Figure [T] For e > 0, r £ <C e, and (m, n) € Z 2 , let 

5ra,n, E = Der E {me, TIE) , 
Qm 7 n,e -^m,n,e ^ ( -F 5 

where D ere (me,ne) is the disc of radius er e centered at (me,ne). Introduce the index set 

I E = {(m,n) G Z 2 : Z? m ,n, £ C u> } , 

and set 

We will assume that the conductivity parameter of Q is of the form 



ri mn\Q s , 



where 

Additionally, we assume 



< k_ < a £ ?rr 2 < k + . (2) 



2tt 2tt 
< — < -£ 2 lnr £ < — . (3) 

7+ 7- 



We consider, for $ b e C 1 ^) and F € L 2 (^), the solution ^ e iT^J?) of 

J -div(a £ V<F e ) = F in i? , 



(4) 



This model, initially introduced by Fenchenko & Khruslov |19j . has been studied in the context 
of homogcnization by several authors [2"1,5,12,9,27J. It is known to have a non-standard behaviour 
when e tends to zero. Namely, the homogenized limit is not of divergence form, and admits a non- 
local term (see Theorem [3] for the precise form of the limit). While it is clear that, thanks to the 
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cllipticity, a global W 12 (f2) bound holds, one can show that except for special boundary data <l>b, 
the solutions X F £ of Q are unbounded in W l( £(Qa) for any p > 2, see Corollary 2\ This makes 
the situation very different from the case of bounded coefficients: Meyers' Theorem 28] shows that 
solutions are bounded in W lo '^ for some p > 2 in that case. In the case of periodic composites with 
bounded coefficients, Li & Vogelius (235J and Li & Nirenberg [21] showed one can in fact obtain 
W l o ° estimates. However, such improved regularity estimates strongly depend on the contrast of 
the coefficients as shown for example in 22,3113]. 

The first goal of this work is to establish interior C 1,a estimates, uniformly in e, for ] F £ away 
from Q £ . We show that in a set f2 £ = lj £ x (—1,1) improved regularity estimates can be obtained. 
The set io £ is "almost" the complement of the high conductivity fibres in the sense that uj £ D £ = 
and wo \ w e tends to zero with e. Introducing V the solution in Hq(Q) of 

-AV = F, 

we show that 

\\# B -V\\ci,* W <C\\F\\ L , {n) , 

see Theorems [T] and [2] for precise statements. 

We can think of several applications for this work. One could for example use this result to 
establish a posteriori error estimates for the numerical solutions of Q . Our initial motivation came 
from a question related to imaging. In a recent work, Ammari et al. [I] showed that the signature of 
small inclusions inside a periodic medium was determined by the effective properties of this medium, 
in the limit when both the period and the size of the inclusions tend to zero. The motivation for this 
study was to provide a mathematical perspective on the so-called "resolution beyond the diffraction 
limit" verified experimentally [23] . 

Recent developments in material sciences, see Bouchitte et al. [7], have shown that very high 
contrast composite materials, with scalings similar to the ones used in this paper, could be used to 
construct meta-materials with particularly interesting properties (such as materials with negative 
optical indices) . Such composites are out of the scope of [1] . This work relies on elliptic regularity 
estimates shown by Li & Nirenberg which do not apply for large contrast. These imaging problems 
involve either the Helmholtz equation, or the full Maxwell equation. We limited our study to the 
case of a real valued conductivity coefficient in this work. Assume a perturbation of small volume 
G £ is located inside fl £ with conductivity £ C°(f2). The conductivity of the defective medium is 

1 in Q \ (Q £ U G £ ) , 
a-e.d = { 71 in G £ , (5) 



For <Pb G C 1 (J7), consider W e ^ the solution of 

/ div(a £ , d VW e , d ) = in n , 
\ W e4 = $ h on 3Q , 

and compare it to the unperturbed solution W e of 

div(a £ VW e ) = in (2 , 
W £ = <P b on dQ . 



(G) 



(7) 



The signature on the boundary of the defect is characterised by the quadratic form 3i e : ff 1 / 2 (9^2) — > 
E defined by 

®e{$b) = [ a £ ®- (W £ . d - W £ ) {s)$ b {s)ds. 
Jon on 
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We show in Theorem [4] that at first order, for G C 1 (.!?), it is of the form 

@ £ ($ b ) = \G £ \ [ MVW* ■ VW,dfi + o{\G £ \), 



n 



where W* is the solution of the homogenized problem (13) associated with This can be seen as 
an extension of [13] to At first order, the signature of the defect is similar to what would be 
observed if defect was introduced in the effective medium, instead of the homogeneous substrate: 
the only difference is the formula of polarisation tensor M . This is what was observed in [5] for finite 
conductivities. 

Theorem [4] has another interpretation, probably of equal if not greater importance for applica- 
tions. Theorem |4j compared with the main result of [1] or [6] shows that impurities in the substrate 
do not affect the overall properties of the composite material more than impurities would affect a 
regular composite, provided these impurities are located in f} £ , that is, not too close to the highly 
conducting fibres but not necessarily at a distance proportional to the size of the microstructure. 
As far as the authors are aware, this is the first result of this nature. If such highly contrasted 
structures are used to manufacture meta-materials, as it is suggested in [7], this result is of practical 
importance. 

Another very related question is the regularity of the solution of a problem similar to Q , where 
cavities replace inclusions. This problem shows a similar non-standard effective behaviour, the cel- 
ebrated "strange term" of Cioranescu and Murat [T7]. This will be the subject of a forthcoming 
paper. 

Our paper is structured as follows. In Section [2] we state our main results concerning the regu- 
larity of |4]) away from the highly conductive fibres. First, we consider for F G L 2 (f2), the solution 
U e G Hl<J2) of 

' -div(a e VU e ) — F in fl , , , 

U E = on df2 . [ ' 

Our result concerning problem Q is Theorem [II We then turn to the boundary value problem 
Q, and obtain an estimate for W e in Theorem - ]^] To highlight the fact that excluding a buffer 
zone around the highly conducting fibres is necessary, for any p > 2 we provide an explosive lower 
bound for the W 1,p norm of W £ in Corollary [2] (see also Remark [3]). This result is a corollary of the 
homogenization result given by Theorem [3] 

Our approach relies on the fact that one can perform a partial Fourier series decomposition of 
U s and W e in X3 . In Section [3] we prove Theorem [l] and Theorem [2] starting from the regularity of 
the Fourier coefficients. Note that translating the interior regularity results from U e to W e is not 
obvious. A natural idea is to study W £ = i]W £ , where 77 is a cut-off function. Provided the cut-off 
function is chosen carefully, this leads to an interior problem in H^{Q) 1 

-div(a £ VW E ) = -2a £ VW £ ■ Vr] - a £ Ar]W £ in Q. 

Let us focus on the first right-hand side term. Remembering that the only bound uniform in e is the 
energy bound, j Q a e \VW e \ 2 < C\\$b\\c l (f))i since Vlf £ is not bounded uniformly in e in any L P (Q) 
for p > 2, we are then left with a problem of the form 

-div(a E VWe) = -a £ f, in H(j(f2) with / a £ f 2 dx < 00. 

Jn 

Unfortunately, without additional information on /, this is not enough to guarantee that W £ is 
bounded in L°°(J7). In Proposition [4j we show that the best one can hope for in this case is 
e\\W £ \\ L ~ <C\\^cr £ f\\ L 2. 
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In Section [4j we show how these results can be used to obtain a representation formula for 
inclusions of small measure located away from the fibres. This part was the initial motivation of this 
work. Because outside of the highly conducting fibres the substrate is homogeneous, we provide a 
self-contained proof. 

We prove the homogenization Theorem [3lin Section [5] and finally turn to more technical results. 
Section [6] is devoted to the proof of Lemma [lj a supremum estimate of de Giorgi-Moser-Nash type 
adapted to our problem. The proof of this Lemma uses a Poincare-Sobolev inequality proved in 
Lemma |3J Section!?] contains the proof of the regularity results for the Fourier coefficients associated 
with U s . Section |8j is the counterpart of Section [jj for W e . Several intermediate technical results are 
proved in Appendix |A"} 

2 Main interior regularity estimates 

The main part of this paper is devoted to the derivation of interior regularity estimates for prob- 
lem Q, outside of the highly conducting rods. 
For a fixed ujq d uj\ (e uj, define by 

cu e T = {x e lui : dist (x, D £ ) > st}. (9) 

Our result concerning problem Q is the following. 

Theorem 1 For n > ; r > Ke 2 ' 1+ ") with r/ S (§ , l), and < v < 2 (n — |), the solution U e of 
^ satisfies 

ll^e- V|| 01 ,„ (ti ,. x( _ £)£ 3 ) < C(n,ri,v)\\F\\ L 2 {n) , 
where V € H^{fi) is the solution of 

-AV = F. 

This result is proved in Section [3] 

Remark 1 It should be noted that, generically, V may not belong to any W 1,P (Q) for p > 6. The 
above result asserts that the difference U e — V enjoys a better regularity in most of the domain since 
\u>i \ cj^| — ;> as e -> 0. 

We use the notation C for various constants in the paper which are always independent of e. 
When appropriate, we highlight the dependence on the parameters appearing in the statements of 
the results by writing C(a,b,. . .). 

As a direct consequence of the above theorem, we have 

Corollary 1 Under the assumptions of Theorem^ we have 

II^IIl-c^x (-lx)) + \\VU e \\l^x(-l,l)) < C{K,rf)\\F\\ LH n). (10) 
Next, we derive an estimate for the solution W e of Q, given by the following proposition. 

Theorem 2 Let Q% = uj% x (—1,1), with I < L. For k > 0, t > ke ! ('+') with -q e (f,l)? and 
d < v < 2 (n — |), the solution W s of Q satisfies 

\\W e \\ i,»(a$ < C(K,r),V,l) \\$b\\ C i(Qy 

This result is proved in Section [3] 

Remark 2 Once again, as e — > 0, we have \u)x \ w tI = therefore the solution of ^ enjoys a 

uniform C ,!/ -bound in almost the whole domain. 



6 



It is natural to ask whether a uniform global C ,!/ -bound exists for J?. To answer this question, 
first we consider the limit homogenized problem corresponding to Q. To state our homogenization 
result, we need to introduce a capacity function. The function c" is defined, for some a > 1, in fi by 



if r = \J [x\ — me) 2 + (x 2 — ne) 2 < er £ 



In r — ln(er e 



ln(e CT /2) - ln(er e ) 

1 elsewhere 



if r e {er e ,e a /2), for (m,n)eJ e , ( n ) 



In fact, we choose a so that e a > er. Note that c £ does not depend on the variable £3, is periodic 
of period [— e/2,e/2] 3 in f2 , and is equal to 1 outside of J? - This capacity function has been 
ubiquitous in the derivation of homogenization results related to conductivities of the form ([I]) since 
its introduction in j!7j . 

We have the following homogenization result, proved in Section [5] Note that a £ is only reinforced 
in the cylinder f2 = cj X (~L, L). 

Theorem 3 Assume in addition to pip, Q), that 

2ir 

lim a E nr 2 = k£ [«-,«+], lim , =76 [7-, 7+]- (12) 

Let <Pb G G (fi). Then, the solution W e converges weakly in to the unique solution W* of the 

coupled system 

-AW m +'y(W m -V,)ln =Oinn 

-k0 2 3 V* + 7 (K - W*) = in !2 = wo x (-L, L) 

W* = & b ondf2 

V*(-,±L) = <P b (-,±L) m w - 

T/ie pair (W*,14) satisfies W* € C ' a (J7)nW^' c p (/?) and 14 € C*' a (A))nCi~ (/2 ) /or some a £ (0, 1) 
and for any p > 2. 

Moreover, the following corrector result holds: 

lim / a e |Viy e -V<(^-V;)-<V^--(l-<)a3V;e3| 2 da; = 0. (14) 
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Because of the non-local nature of the problem it solves, it is not clear that W* is analytic inside 
Qq. The following proposition shows that it enjoys partial analyticity, that is, with respect to the 
X3 variable. 

Proposition 1 Under the hypothesis of Theorem^ the solutions W* and V* of (13 1 are analytic 
with respect to the X3 variable in Qq. 



This result is proved in Appendix |A.1| We are now in position to answer the global regularity 
question. Theorem [3] and Proposition [1] imply that one cannot expect an s independent bound for 
the sequence W e in a space better than H 1 (f2) in general, as the following Corollary shows. 

Corollary 2 Let <L> b € C 2 (S2) be such that the sets 

S±($b) = W € ujq such that A 2 @b(-, ±L) = 0} , where A 2 = + <9f 2 , 

have an intersection of zero measure. Then, the sequence W £ is unbounded in W^(Qq) for any 
p > 2. More precisely, for any non empty open set fl' C £2q, there exists a positive constant C(fi' ', <Pb) 
independent of e such that for any p £ (0, 1) and p > 2, 

liminf (^-^IIVVMWaQ,)) ^ i 1 P) C ^ ^b) . 
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Proof Working with a subsequence if necessary, we can assume that (12) holds. 
Fix p G (0, 1) and let O e be the open subset of S7o defined by 

O, = {x e u : < dist (x, D B ) < er§} x (-L, L), (15) 



By the definition ( 11 ) of c° combined with ( 12 ), and by the periodicity of c% and lo e in f2o, we have 

lo. |VcJ| 2 - lim - / - U 7(1 " P) weakly- * in A1(/2 ), 

e->o y e 2 J £Tf r In (2e 1 " CT r e ) J 

which by virtue of the continuity of V* — W* implies that 

la. |Vc^| 2 (K -W,) 2 dx 7(1-/3)(K -W*) 2 weakly-* in Af(f2 ). (16) 

Let us show that this last term does not cancel on any non-empty open subset Q' C J?o- By 
contradiction, suppose that the continuous function V* — W* = on a ball Br(p°), centred in 
P° = (PiiP^iPs) an( l of radius R > 0, such that Br(po) C X2o- This implies that for example that for 
any x' G C fl , with C R = (pj- R/2,p\+R/2) X (pg — iJ/2, p§ + J2/2) , the function H^fx', •) -V*(a/, •) 
is identically zero on (pi - y/ZR/2,p% + V3R/2). As we have shown in Proposition 1 this function 
is analytic on (— L, L) therefore — V* on Cr x (— L,L). The system (13) then shows that 
9^V* — AV* — AW* = in Cr x (— L,L). Hence, taking into account the boundary conditions of 
( 13 ) we would have 

W,{x) = V,{x) = ^p-$ b (x',L) + ^-^$ b (x',-L) for any x G Cr x (-L,L). 

Since W 7 * is harmonic in Cr x (— L, L), it follows that both functions <£&(•, ±L) are harmonic in Cr, 
therefore \S+($b) H 5_(^b)| > |Cr| > which contradicts the hypothesis. 
Next, since VW,, d 3 Vi G i 2 (/?o) and < c% < 1 in /2 , w e have 

lim / |<V^ + (l-c^)93Ke 3 | 2 da; = 0. 
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Hence, by the convergence (14) and the fact that a e = 1 in O e , we deduce that 
lo J VW e ~ V<£ (W» -V,)\ 2 — > strongly in Lj, c (ft„). 



This combined with convergence ( |16[ ) implies that 

lo e |VW e | 2 -» 7(1 — /o)(F* — H 7 *) 2 weakly-* in M(Q ). (17) 

To proceed, fix Q" @ fi' C fio and select a smooth cut-off function £ G C^(fi') such that 
< C < 1 and C = 1 in i?". Note that, for some C x > 1 

CfVf" < |O e | < Cir 2p . 



Thus, (17) shows 



7(1 -p)(V, - wg 2 dx < / 7 (l-p)(K- Wj'Cdx 

IQ 



= lim / \VW e \ 2 Cdx 



>0 



< liminf / \VW F \ 2 dx 



o e nn' 



< liminf (|O e r- 2/p ||VW e ||L (fi , We) ) 

< d hminf (^-^llVVMLtrnQ.) 



<s 



Since this is valid for all O" <s f2', we conclude the proof. □ 

Remark 3 Note that Corollary[2]shows that gradient blow-up occurs outside of the highly conducting 
rods. An variant of the proof shows that one can obtain the following estimate 

liminf (r e 2Kl - 2/p) ||VW e ||^ (fi/ ^ r) ) > (p'-p)C(Q',$ b ), 
where < p < p' < 1 and 

QP' = {x e uj : dist (x, D e ) < er( } x (-L, L), 
which shows that the blow-up is not localized on the surface of the rods. 



3 Regularity estimates in two dimensions for weighted equations and 
three-dimensional consequences 

In this section, we state the regularity results we have obtained for two-dimensional companion 
problems of problem Q . To highlight the sharpness of the supremum estimate provided by Lemmajlj 
which provides an L°° upper bound of the form C{s)e~ s for any 2 > s > 1, we construct in 
Proposition |4] an example where the L°° norm is bounded from below by CeT 1 . Then, we prove 
Theorem Q] and Theorem [2j 



3.1 Towards the proof of Theorem [T] 

Let us now turn to the proof of Theorem [T] Because none of the coefficients depend on X3, we can, 
as in jTU] , reduce the study of this three-dimensional problem to that of a two-dimensional problem 
by separating variables. We write 



00 

U e {x',x 3 ) - ^2u £in (x') sin (Jy(°T + ^) ' 



n=l 
00 



F(x',x 3 ) - ^2f n {x') sin 



717T ,£3 
~2~^ ~L 



(18) 
(19) 



and ^ becomes 



-div 2 (a e V 2 u e 
u E .„ = 



4^2 Oe W e>; 



fn in LU , 

on duj . 



(20) 



In the above, div 2 and V 2 denotes the horizontal divergence and gradient operators. We are thus 

2 

led to consider, for A > f^j, 



— div 2 (a e V 2 u £ ) + A a e u £ = f + a e g + div 2 (h) in uj , 
u e = on du 



(21) 



where / e L 2 (uj), ^/a^g £ L 2 (uj) and h 6 L°°(lu) 2 . Problem (21) is equivalent to Problem (20) if 
g = 0,h = 0; the additional a e g + div 2 (h) term will prove useful for the study of boundary value 
problem (17 1. 

Section 7] is devoted to the proof of the following result. 
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Proposition 2 Assume that A > Ao > 0. For n > 0, 77 6 (0, 1), and t > ke^+i) i/ie solution u E to 
(21 1 enjoys the following bound 



< 



C(A ,/t, 77) 
A? 



A* 



a (ll/IU 2 M + HV^SIU 2 ^) + ll h llL-( w )0 + |K|U«M 



where v G Hq(lj) is the solution of 



Furthermore, for < 1/ < 77, 

C(A ,k,?7,^) 



-/Au + At) = / + .g + div 2 (h) 



< 



A2 2 



A 



a (ll/IU 2 (w) + WVa^ 9\\l 2 (uj) + WHl^(u) 2 ) + 



Throughout the paper [-Jo* will be used to denote the ^-Holder semi- norm. The second ingredient 
for the proof is a supremum estimate for u e . 

Lemma 1 Assume that conditions (|Tj) , ([2j) and Q hold. Assume that (p e 6 Hq (lj) satisfies, 

- div (a e Vtp € ) + \a e ip e = a e / e + div (a £ h) . (22) 
Then, for any l<a<2, 0</3<l — a/2, we have 

.. C(a,(3) /.. ., 2+1 11 ,11 \ 

This variation on the standard de Giorgi-Moser-Nash estimates is proved in Section [6j As a direct 
consequence we have the following proposition. 



Proposition 3 Assume that A > Aq > 0, g = and h = 0. We have 



n n ^ C(X , q) 



(23) 



for any q > 2. 



Note that Lemma [I] does not show that H^Hl^u) is bounded independently of e. The following 
counter-example documents the sharpness of our estimate, as it shows that it must be at least 

Proposition 4 Assume that conditions (JT1), |2| and Q hold. Let ip e £ Hq(ui) be the solution of 
(22 1. Let lrj g e be the indicator function of the disk -Do,o,e C R 2 centred at the origin and of radius 
er e . Set f £ = e~ 1 1d e and h = 0. Then, we have 



> 



1 K- 



where lim e _j. <?(1) = 0. 



(24) 
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Proof of Theorem^ Fix some I € (0, L) and let Q% — uj^ x (— Z, Z). It suffices to provide the required 
estimates in fl € T with a constant that is independent of I. 



Together with ( |18[ ) and ( 19 ) we consider the Fourier decomposition 

oo 

'717T ,X3 



n=l 



Note that £/ e — V is harmonic, thus regular, in Q \ Q e . Therefore, by Dirichlet's theorem on the 
convergence of Fourier series, the Fourier expansions of U £ — V and V(Z7 £ — V) converge pointwise 
to U e — V and V(U e — V) in fl%. Proposition [2] and Proposition |3] combined with the choice q = | 
show that 



C(/c 7]) 

IK,n - Ur»|U°°(u<?) + ll^U^n - V 2 « n |U°°(u>?) < 2 ^ ll/n|U 2 (o.). 

C(k 77) 



with < v < 77 < 1. For 7/ > 4, (25) shows that the sums 



h(x' ,Xz) = ^2 ( u e,n( x/ ) - Vn{x')) sin 



717T /X3 



1 , 



n7r ,£3 



^2(^,^3) = ^2 (V2U e , n (x') - V 2 w„(a;')) sin 

71=1 

are absolutely convergent, 

00 ^ 

||A||i«(fJf) + ||-f2|U«(fi ? ) < C z2-Zpj\\fn\\L'(w) 

71=1 

<cii^iu^) (E^) 2 <Aii^i 

\n=l / ' 4 



L 2 (J?)- 



This shows that 



\\Ue ~ V\\ Lae{nm) + ||V 2 E4 - V 2 ^|| ioo(n;) < C\\F\\ LHn) . 



We turn to the Holder estimates of V2U £ — V2V By (26) with 77 > | + | we have, 

00 

l/sCa;',^) - h{y',x 3 )\ < J2 |(V a «e,n(a!') - V 2 u„(a:')) - (V 2 u e ,„(y') - V 2 «„(i/'))| 

n=l 
00 

71=1 

< — £-^Mi*(fl)l^-!/T. 

f- 2 - I 



(25) 
(26) 



11 



By (25 1 but with r\ > | + \ we obtain, 



UTT , X 3 

1 2 < I. ' ' 



IW.a*)- W,y 3 )l 

OO 

< l v 2M e ,«(y') - V 2 f„(y')l 

oo 1 

- C 12^j H/«lli 2 l na; 3 - ny 3 \" 

n=l 

< _°_ 1 \\F\\ L2{n) \x s -y 3 \". 

From these last two estimates, we deduce 

[V 2 U e -V 2 V] c „ {acT) <C\\F\\ L * {Q) . 



Next, we estimate d X3 U e — d X3 V. Choosing r\ > § in (251 shows that the sum 



n7r , x 3 



7 3 = ^n (u e , n (ar') - v n (x'))cos } , 
converges absolutely with bound 

oo 1 c 

\\h\\L°°(Q*) < r^Fill/™IU 2 (") ^ r3ill F IU 2 («)' 

n=l ' 4 

This proves that 

\\d X3 U E -d X3 V\\ LOO(neT) <C\\F\\ L2(n) . 



For the Holder estimate, we proceed as before. First, with rj > 4 in (25) we obtain 



\h(x',x 3 ) - I 3 (y',x 3 )\ < ^2n\(u e>n (x') - v n (x')) - (u e>n (y') - v n {y'))\ 

n=l 

oo 



< 



n=l 

-^3||F|U 2(fi) |^-y'|. 



Also by ([25} but with 77 > § + |, 



IW.^s) -h(y', 2/3)1 





^3 


n)) 


cos(- 











V 2 v i 



00 1 

\\ F \\L2(n) \%3 - y 3 \ v - 



< 



n=l 
C 



n - - - - 

''2 4 
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We have shown 



[d X3 U e -d X3 V] c „ m) <C\\F\\ v . 



(«)' 



which concludes the proof. □ 



Proof of Proposition The bound on f E in (24 1 is straightforward, from the definition of a E , ([!]), 
and ([2]). Integrating (22) against ip E and using ([2J), we find 

/ a e |V( / 9 e | 2 da;+ / \a E ip 2 E dx = I a E f E ip E dx = / a E (p E e^dx < k + e \\ip E \\ L00 ^y 

Let Y = ( — 1/2, 1/2) 2 be the unit period cell, and set 

9e = e- 1 (l-c 1 s )l eY eH^u J ), (27) 
where c\ is defined by (11) for c = 1. An easy computation using ^ and ([3]), shows that 

/ a E \Vg E \ 2 dx + / \a E g 2 dx = / | Vc* | 2 e~ 2 dx + / A a E (1 - c\) 2 e~ 2 dx 

Joj Juj JsY JeY 

< 7+ + Ak+ + o(l). 

Then, the Cauchy-Schwarz inequality combined with the two previous estimates yields 



a e V^ e -Vg e dx + I a £ Lp e g e dx<[K + e\\Lp e \\ Lao{ul) y(^+ + \n + +o(l)) 



On the other hand, integrating (22) against g E: we find 



a E V(p E ■ X7g E dx + / a E (p E g E dx= / a E f E g E dx— / a E e 2 dx > k_ 
Combining these two inequalities provides the announced bound. □ 



3.2 Towards the proof of Theorem [2] 

Let us now turn to the proof of Theorem [2j Our approach is to decompose the boundary data #b 
into three parts. We first define 



<>l(;Xs) = g (M; L) - -L)) + \ (M; L) + -L)), 



(28) 



which agrees with <?;,(•, ±L) onwx {±L}. We now freeze the variations of 4>l in the a/-direction on 
the highly conducting fibres by decomposing <f>L into 



where 



•)i=(t>LC a E + ^2 (l-c E )l {l ( Xl . X2) -( m , n)El < E /2}(x)(j) L (ms,nE,x 3 ), (29) 

(m,n)e/ e 



the capacity function c E being defined in (11). We finally define </>o to be the trace of on the 
lateral boundary of the cylinder SI, 



on dfi. 



Note that 4>q vanishes onwx {±£}- 

The properties of <f)\ and 4>2 we will use are described by the following proposition. 
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Proposition 5 The function (f>i given by ( 29 ) 

— belongs to C°°([-L, L]; C ' 1 ^)), 

— is independent of Xi and X2 inside each connected component of Q £ , 

— satisfies d^(pi = in fi, 

— is globally Lipschitz, 

\\<t>i\\c°,i(n) ^ C\\(j} L \\ cl ^y 
The function </> 2 — <Pl ~ 4*1 where 4>l is defined in ( 28 1 

— is globally Lipschitz, 

\\fo\\ c o,i(n) ^ CMdc^ti)' 

— is supported by 

— s x £ Slo s.t. min \(xi,X2) ~ (me,ne)\ < e a \ . 

\ (m,n)el c ' J 

— and satisfies 

\\<hh-(0) < Ce^U^^y 
The proof is straightforward and given in Appendix |A.2| 

Proof of Theorem^ We will study three boundary problems separately. For i = 0, 1, 2, we introduce 
W Ey i, the solution of 

' div (a s VW Bii ) = in Q 
W e ^i — cf>i on df2 

By linearity W £ = W efi + W eA + W e>2 . 

Consider W e fl — Co W Ey o where Co is a cut-off function Co € C^°(w3), such that 



(30) 



< Co < 1 and Co = 1 in w 2 , (31) 

with W! g w 2 6 W3 g lo. Since Co = 1 orL w Qj W e G Hq([2) satisfies 

- div (0= W e , ) - -2VCo • W £)0 - (^Co) W £ , in fl. (32) 

The maximum principle shows that || W e ,o || L°°(n) < Il0o||c°(r2)- Furthermore, on = (cj \ cjo) x 
(— £, L), u e = W e fi is the solution of 

Zlw £ = in 
u e = on (w \ ujq) x {-£, L}, 
u e = W £fl on <9(w \ uj ) x (— L, L). 

Classical regularity estimates then show that 

H^oli c2( ( ^ 2)xW) ) <C{u 2 ,U 3 ,l)\\WeA L oo {Q y (33) 

Thus, the right-hand side of ([32]) is in C 1 (/5), and we can apply Theorem [l] to conclude that for 
t > Ke^a+^y w ith 77 e (§, l), and < v < 2 (77 - §) , we have 

ll^.o -^|| O x., {n . ) <C||W e , || ioo(n)> 

where Vo € Hq(Q) is the solution of 

-AV = -2VCo • Vlf £ , - (4Co) VF £ , in 12. 
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Since the right-hand side is C 1 (/2) , Vq enjoys interior regularity, and we have obtained 

\\W E , \\ c ^ m) <C||<M c o (fl) . (34) 
Next, we consider W e> \. We note that for any v £ Hq(Q), 



a s V (W £) i - fa) ■ Vu = - / a £ Vfa ■ Vv, 

Q 

(1 - a e ) V^i • Vu - / V^i-Vv, 
(1 - o e ) <9 3 </>i • <9 3 u + / Afav, 



= / A£iu. 
J n 

We used first that 4>\ is independent of £1 and x 2 in each connected component of Q £ , and then that 
it satisfies d 33 4>i = (see Proposition [H]) . In other words, introducing W St i — W e ,i — <pi, we have 

- div (o 6 W Ei i) = Afa £ Q. (35) 

Following the strategy used for the interior source problem, we thus introduce V\ £ Hq(Q), the 
solution of 

-AVi = A<f>i in Q, 

In Section [sj we will show in Proposition [9] that, for r > ke'C+'I , | < 77 < 1, < v < 2(j] — ~) that 

11^,1 - ^lllci^fl.) <C||#6|| C i(fl)- 

On the other hand, as V\ + 4>\ is harmonic in f2, we have 

As W^i = (We,i — Vi) + (Vi + </>i), we have thus shown that 

Finally, for W E 2) we shall consider instead the function W Ey2 = C where £(2;) = p(x3)Co and 
£0 is the cut-off function given by (31 ), and p £ C%°(—1', I') with I < /' < L is such that 

<p < 1 andp= 1 in (-1,1). (37) 

We compute 

- div (a e VW e . 2 ) = a e (2d 3 (d 3 W s , 2 + d 33 £ W s , 2 ) + g e , (38) 
where g e has support in (0/3 \ 0J2) x ( — I', I'), and is given by 

g e = (d n C + d 22 () W e , 2 + 2d 1 Cd 1 W efl + 2d 2 (d 2 W e , 2 . 
By the standard maximum principle, 

l|W 6 , a || LOO(n) <||^|| £ . (n) . (39) 
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Just like W Sj o enjoys additional regularity away from the support of the rods as shown by (33 1, 
we have, for any I' < I" < L, 



(40) 



Next, notice that d X:j W £> 2 and W £t 2 are both a 6 -harmonic, a standard energy estimate shows that, 
for any S2' <s Q" <s ,!?, we have 



/ a e \V(d 3 W e , 2 )\ 2 dx<C [ a s \d 3 W s , 2 \ 2 dx<C [ a £ \W £ . 2 \ 2 dx < C \\W ea \\ L , 
J n> Jn" Jn 



■(«)• 



Note that the regularity estimate ( 40 ) yields in particular that 

belief) < C\\<h\\L^{n) ■ 



The companion problem to ( 38 ) is 



= 2 ^Ca 3 W £>2 + d 33 (W £ , 2 + g e in Q , 



(41) 



(42) 



(43) 



Thanks to (40 1, the right-hand side of (43) is in C 1 (l?) uniformly in e, and therefore V e2 enjoys 
interior regularity. Proposition 10 proved in Section 8 on the regularity of the difference W £i2 — V £ . 2 

1—7} ' 

implies a regularity result of W £ . 2 . For n > 0, r > k£ ! ('+') with rj G (|, l), and < v < 3 (77 — |j, 
the solution W e ,2 of (30) satisfies 



<Ce°\\$ b \\ c , {ny 



(44) 



The claim of the Theorem is a consequence of the three bounds (34), (36) and (44). □ 



4 Application to structures with defects of small volume 



We now consider the case when a defect of small volume is present in the medium. In that case, the 
conductivity of the defective medium is given by ((5|. We assume that the defect of support G £ stays 
away from the high conductivity fibres. To fix ideas, given < I < L and uj± such that luq eg ui <g w, 
suppose 

G £ C f2 £ = {x e ^ : dist (x, D £ ) > e 17/le } x (-1, 1). (45) 

Note that the set fl £ grows as e tends to zero. Furthermore, in the language of Section [2] Q £ = fl e T 
for r = e 1 / 16 , k = 1, rj = 7/9. This guarantees a C 1,u estimate with < v < 1/18. 
Throughout the section, we will assume that ( 12 ) holds, namely 



lim a F 7rrv 
£^0 £ 



K £ [K— , 



lim 

e^O e 



2tt 
2 l mr 



7 e [7-.7+J- 



When e — > 0, Theorem |3| sh ows that the sequence W £ converges weakly in Hq(Q) to the solution 
W„ of the coupled system Jl3\ when <P b e C 1 ^). 

We proceed to derive an asymptotic formula for the difference response operator (the difference 
of Dirichlet-to-Neumann operators on dSl) 

H 1/2 {dQ) -> H- 1/2 (dQ) 
d 



$h a. 



dn 



(W £td - W £ ) 



on 
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As this operator can be recovered by polarisation, we limit ourselves to the study of the quadratic 
form 

& e : H 1/2 {dQ) -> R 

<P b ^ [ a^{W E . d -W £ )(s)Ms)ds. (46) 
J do dn 

We can now state the main result of this section. 

Theorem 4 Assume that \G S \ — » 0, and \G £ \~ 1 lG e converges weakly-* to a Radon measure /i on Q. 
There exist a subsequence, still denoted by e, and a matrix-valued {unction M £ L 2 (J?,^) 3x3 such 
that the bilinear response form £$ £ given by ( 46 1 has the following asymptotic form 



M e {$ b ) = \G £ \ / AfVW 7 * • VW* dfi + o(\G e \) for any$ b €C L (n), (47) 



n 



where W* denotes the solution to the homogenized problem (13) with boundary condition <P b , and 

lim sup = 0. 

In addition, the matrix-valued function M is symmetric, independent of<P b , and satisfies 
( 7l -l)minfl,~) <M(y)£-£< ( 7l - 1) max ^1, ^ fi-a.e. 



Remark 4 Due to the definition (451 of fi £l th e measure /x has actually support in uj\ x [—1,1] C= i?. 
This gives a sense to the integral term of (47) since € Cl oc {Q). 

Remark 5 We have characterised the signature of a defect provided that the boundary condition <P b 
is in C 1 (^7). The natural space for <2> b is H^ 2 (dn), but we do not know if our result holds in this 
space as well. 



Our strategy is inspired by (131 11] . The first key element is the following Lemma. 
Lemma 2 There exist a positive constant C independent of <!> b and e such that 

||V(W £ - W e , d )\\ L 2 {n) + [G^^WW, - W e , d \\ L 2 [Q) < CIGJ 1 / 2 H^Ho^fi). 
Proof Write X e,d = W S - W e . d G H^(Q), we have 

Aw{a e4 V X e.d) = div((a £ , d - a e )VW £ ) = div(( 7l - 1) l Ge VW E ) 
Thus, by testing against Xe,d, we get 

\WXe,d\\L-(n) < C\\VW s \\ LHGe) < C \G e \ x / 2 \\VW e \\ L <*> [Ge) . 
Next note that, thanks to Theorem [2j we have 

||VWe|U~( Gl ) < \\W e \\ c ^(n e ) < C\\<P b \\ (48) 

with v = 1/20. Thus 

l|v XM || L2( , 2) ^cig^w^w 
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This proves the first half of the desired estimate, namely 

||V(W e - W e>d )\\ LH a) < C |G £ r/ 2 W$b\\cHny 
To prove the second part of the estimate, introduce the solution F e e Hq(S7) of 

-div {a e WF e ) = Xe,d ■ 

Thanks to Corollary [l] we have 

\\VF e \\ L6{ae) <C\\x s , d \\L*(n). 

On the other hand, note that 

div(a e VxM) - div((a e . d - a £ )VW e4 ) = div(( 7l - 1) l Ge VW e , d ). 
After an integration by part, we see that 

IIXMlli^fl) = J a e VF e ■ Vxe,ddx 

(71 - 1) l Ge VF £ • VW £ , d dx 

(71 - 1) VF e ■ Vxe.d dx + f ( 7l - 1) VF £ ■ WW e dx 

<C||VF £ || W) (||VxMll i6 /5 (Ge) + ||V^ e ||^ (Ge) |G e | 5 / 6 ) 
< C\\xeAL* { n) (llVxe.dllL^^IG,! 5 / 6 - 1 / 2 + ||VW e |U~( Ge )|G E | 5 / 6 ) 
<C\\xeALHn)\G e \^\\<P b \\^ {n) , 



G, 



where we have used (49 1 in the penultimate inequality. The proof is complete. □ 



The second ingredient is a pointwise uniform estimate on W E — W*. We use the following extension 
of Theorem[3]with varying boundary data, the proof of which consists of a straightforward adaptation 
of the first step of the proof of Theorem [3j 

Theorem 5 Assume that the conditions and |7^) hold. Consider a sequence ip e 6 G 1 (i7) 

which converges strongly in C°(f2) to some if € C 0,1 (f2). Then, the solution W e 0} problem with 
the boundary condition = (p e converges weakly in H 1 {Q) to the solution W* of ( |13| with the 
boundary condition = cp. Furthermore, W* € W^(f2) and V* € ^^(fio) f or any p > 2. 

Proposition 6 The solutions W E and W* to respectively ^ and ( 13 ) satisfy the following conver- 
gence in the set Q e of (45), 



lim sup ||VW e - VM / *|| G o,i/2 0(fie) = 0. (50) 



$ b e C\f2) 
Wb\\c H n) < 1 
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Proof By contradiction, suppose that there exist a constant Co > 0, a sequence (e n )neN, and a 
sequence ip n £ C 1 (/2) such that 

||¥>nll i(fl) - 1 > 

and 



|VW„-VW*,„| 



> Co, 



where W n is the solution of 

div(a e „VW n ) = in J?, and W n = y>„ on dfl, 

and W*n is the solution of 

- zAW*,„ + 7 (W*,„ - 14,„) lfl,, = in Q 
- Kd$ 3 V* t n + l(V*, n - W*,n) = in J? = w x (-L,L) 

W* <n — (fin on 9/2 
K,n(-,±£) = ^n(-,±i) in uj . 

The regularity Theorem [2] shows that there exists a constant K, independent of n, such that 



ivw, 



n|| C o,l/l9( fie ) 



< K \\<p n \\ cl{ n) < K 



Then, by virtue of the extension theorem for Holder functions (see, e.g., |30j ) there exists an extension 
£„ of VW n to the set Q\ = ui x (— Z, Z), such that 

||Cn||c°. 1 /i9(f2 1 )3 < if, 

with the same constant K. Note that the embedding C 0,1 / 19 (J?i) C 0,1 / 20 (J2i) is compact. There- 
fore we can extract a convergent subsequence £ p (and the associated sequence tp p ) such that 



0. 



Thanks to the Ascoli-Arzela Theorem, we can extract yet another subsequence £ g , tp q such that 

II — vllcofiJ) — > f° r some V G C 0,1 (f2). 
Now, the homogenization Theorem [3] shows that 

VWg ->> VW* in L 2 (^) 3 , 

where W* is uniquely defined by 

( - AW, +7(W» - 14) l flo = in 12 

- k9| 3 K +7 (14 - W*) = in 0) = wo x (-£, L) 
W* = <p on 9i7 
T4(-,±L) = ip(-,±L) in wo- 

A uniqueness argument shows that VW* = £ on In particular, 



\VW q -VW4 co , 1/20(Qe) = o(l). 
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On the other hand, thanks to Theorem [3j the linearity of the homogenized system W* jq — £ 
W 2 ' A (Qi) and the Sobolev embedding Theorem we have, 

||VW% - VM / *|| C 0,l/20 (f2e) < C\\<p q - <p\\ C o( Q) = o(l) 

for a constant C independent of q, ip and e. We therefore have obtained that 

C < \\VW q - VW*,g|| CO , 1/2 (fiE) = O(l), 

which is a contradiction. □ 

The third ingredient is an asymptotic representation formula for V (W £y d — W £ ) in the set G £ , 

Proposition 7 (a) Let v\ andv\ d be solutions to Q and pj, respectively for <5>\ = Xi. Then, there 
exist a subsequence, still denoted by e, and a matrix-valued function M* £ L {Q, /i) 3x3 , such that 

for all ip € C°(i7). Furthermore, the matrix M* is symmetric, and satisfies 

(71 - l) 2 

< M*(y) < ^ '- fx-a.e. (52) 

7i 



(b) We have, for all w € W 2 ' i {f2), 



(a £ - a M ) V (W e , d - W e ) -Vwdx = \G £ \ J M*VW, -Vwd[i + o{\G £ \) \\w\\ w2A(n) , 
where o(\G £ \) / \G £ \ converges to zero uniformly for \\@b\\ci(a) — 1- 

Remark 6 It is natural to ask if the extraction of a subsequence is necessary, even in the case of 
bounded conductivities. The answer is positive. In the case of domains G £ shrinking to a point, there 
is an entire set of possible polarisation tensors for a given constant 71, delimited by the so-called 
Hashin-Strikman bounds, see e.g. [SS p g i lTi P ] . 



Remark 7 We can use a compactness argument to derive from (51 ) that 



lim 

e->0 



sup 

IWIc«-(fi)<l 



P. 



\G F 



dx; 



ip dx 



= 



for any < v < 1. We remind the reader of Remark [4} [i has compact support in fl. Indeed, if this 
fails for some v, then there exist some Cq > and a sequence ip n € C°^(f2) such that (along some 
subsequence e —> 0) 

ig d (Kd- v i) f 

ip n dx - I M*. ip n dfj, 



\G £ 



dx 3 



> c . 



On the other hand, as the embedding C a '"(f2) C°(f2) is compact, we can assume that ip n — > V* 
in C°(i7), which implies a violation to (51). 
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Proof The proof of this proposition is a variant of that of the main theorem of [T3] ■ We provide the 
proof here for sake of completeness, adapted to the case at hand. 

(a) Introducing v\ and v\ d solutions of solutions to Q and (jfi]), respectively for a given <l> l b — Xi, we 
know thanks to Theorem |3] and Proposition [6] that 



lim 



\H - vl\\ L 2 (n) + \\W E - V«* 11(70.1/30(0.) 



0, 



where v\ is the solution of W* of (13) corresponding to #b = x^. It is immediate to verify that 
v\ = Xi. In particular, we have 



\W 1 -^(Q e ) 



< C Vv 



\c°' 1 / w (n e ) 



) =o(l), 



(53) 



Here and in the sequel, o(l) denotes any quantity going to zero with e, independently of <&b- On the 
other hand, using Lemma [2] we have 



< C \G E \ 1/2 and lit;* 



Applying the Cauchy-Schwarz inequality yields 



,5/6 



1 

|Gl 



G c 



|V< rf -V<| < j±- \G t 



1 1/2 



v< d -v<|| L2(Ge) <c, 



thanks to (54). We may therefore extract a subsequence, still denoted by e, such that 

IGep 1 1 G , (1 - 7l ) A « rf - V\) dx — ► dMtf, 



(54) 



(55) 



(56) 



where -Mij, 1 < i, j < 2 is a Borel measure with support in uq x [— 1,1]. The above convergence 
results hold in the weak-* topology of M. (ui x [— I, I]) (the set of Radon measures on Q x [ — ^ (])• 
Note that 1 - 71 € G°(/2). We see that, for any / € G°(J?), 



lim IGJ" 1 



(l-7i) 



G e 



d 



d 



dXj M dXj 



v> fdx 



s v II 1 ~ 7i||i°°(«) M _ j 
< hmsup - |VUg 



li 2 (G e ) 



16* 



fdx 



e\JG 



1/2 



Since |G £ 



fdx 



/ d/x, we conclude, using (54) that 



< G 



l/| 2 ^ 



«0 



1/2 



(57) 



It follows that dAiij is absolutely continuous with respect to /j, and thus, there exists a 3 x 3 matrix 
valued function M* g L 2 (f2, fj,) 3x3 such that 



fdMi 



MZJdfi. 



O 



Let us now turn to the properties of the matrix valued function M* , following |13j . The matrix M* 
is characterised by 



1 

a 



e\ JG 



(l- 7l )V« d - vD-Vx^dx + oil) 
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for any ip € C° (Q) . Note that v\ d — v\ is the solution in Hq{Q) of 

div (a M V - vi)) = div ((1 - 7i)lG.V«j) ■ (58) 
Using this identity together with Proposition [6] and Lemma [2j we obtain 

Mfjipdfi = -L / (1 - 71) V (« d - V) • V^dx + o(l) 



I G £ | J G 



^ aE d V (« d - vl) if>) ■ V (< d - vi) dx + o(l) 
J- / a e , d V« d -^)-v(< d -^)^^ + o(l). 



|G £ | j fi 



Under this last form, it is apparent that M 4 * = M* { . Furthermore, given £ € K 3 , introducing 
V e = Ya=i(K d - v l)&i tnis last identity yields 



[ M*£ • £^ = I Y\ Mijti tiMli = 7^1 [ a £ . d | W e | 2 V da; - 

Jfi 1< - - <3 l^el JQ 



o(l). (59) 



This shows that M* > 0, //-almost everywhere. Alternatively, from (581 we derive, using Proposi- 
tion [6] and Lemma [2] that for ip > 0, 

a M |VU £ | 2 V^= / a eA VV e -V{V^)dx + o{\G £ \) 



jT ((1-7i)1g.) fe Vu ^l -V(U^)d.x + (|G E 



((1 - 7i)lf?J^ ' W e da: + o(|G e |) (60) 

< (^fJ^—^l G M 2 ^dx\ 1 a M |W e | 2 ^dz) ' +o{\G e \), 

where we have used we have used rt53b to derive Mm. This shows that M* < ( 7lj , //-almost 
everywhere. 

(b) Let us now show that 

IGef 1 / (a E - a e4 ) V (W e ,d - W e ) ■ Vwdx = |G S | -1 f M*VW, ■ Vwd^i + err, (61) 
Jn Jn 

where err satisfies 

|err| = o(l) \G e \ \\^b\\c l (n)\\ w \\w^W 



Identity (61 1 can be interpreted as a 'separation of scales' result. On the left-hand side, the per- 
turbation induced by the defect is present both in the 'microscopic' term a e — a ed and in the 
'macroscopic' term VW £i d — VW e . On the right-hand side, the 'macroscopic' term VIU* is inde- 
pendent of the defect. To prove this result, we follow the adaptation of the oscillating test function 
method in homogenization of Murat & Tartar [3T] proposed in [T3j . In the following computation, 
v l d ~~ v l pl avs the role of the corrector function in homogenization. By a chain of integration by 
parts, we will transfer the derivatives from V(W £j( j — W £ ) to V(w* d — v\), and then pass to the limit. 
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To express the left-hand side of (61 ) in terms of M* , we introduce v\ d — v\ in the computation 
as follows. Using Einstein summation convention for the index i and noting that v\ — Xi, we have 

f ■ dw 

(a e -a E , d )V(W e , d -W e )-Vwdx = / (1 - 7x ) V (W e , d - W s ) ■ V< — dx (62) 

(1 - 7l ) V (W eA - W £ ) ■ Vv\ -^dx + em, 
with 

|ern| < C||V«* - V<||co W ||V™|U~(G e ) |G £ | 1/2 ||V {W s , d - W e )\\„ {Gt) 
= \G e \ \\$b\\ci{n)\\w\\w^(n), 
thanks to ( 54 ) and Lemma [2j Continuing the transformation, we write 

f (1 - 71) V (W M - W e ) ■ Vvl ^dx = [ (1 - 71) V ( |^ (W £ , d - We)") • Vujdx + err 2 (63) 

a M V (W E , d - W e ) J ■ V « d - vt) dx + err 2 , 



with 



|err 2 | <C\\\7vl\\ c o {n J\7 2 w\\ L 2 {Gs) \\(W £ , d - W e )|| £3(Ge) <C|G e |" ||#6||oa (fl) IMI^(fi) 

= o(l) |G £ | ||^b||ci(^)IMIw 2 > 4 (fi)> 
thanks to Lemma [2] and Holder's inequality. To remove the term, we write 

° mV {We4 ~ We) ) ' V ( V *> d ~ ^ dX (64) 

a e ,dV (W M - W E ) • V (< d - w«) J dx + err 3 

(o e - a M ) VW e • V f — - «*) J + err 3 , 



with 



|err 3 | <C J \V 2 w\ ( \W s , d -W £ \\V « d - + |V (W M - W e )| |< d - «*| ) ds 



< C||^ll^,4 (r2) |G £ | 12 ||*6||ci(fl) 
= o(l) |G £ | ||<Mc?i(n)IMIw*.*(rt), 
using Lemma [2] together with the interpolation inequality 

ll/IU*(«) < ^II/IIl^II/IIl^xj)- 
Expanding this last expression, we have 



(a e - a e . d ) VW E • V (J^- « d - dx (65) 
(1 - 7i) VW e ■ V (vi <d - v l E ) ^dx + crr 4 , 



with 

1 1/2 



|err 4 | < C||< d - «'|| £a(G .)|| VW e \\ L ~ {G .) |G e | i/z || W || 
= o(l) |G £ | ||<£ 6 || c i ( ^)|Mlw^(r2)- 



Using the convergence of W e to W„ , we get 



fill) 

^ (l-7i)VW e -V« d -^)— dz 
1 VW, • (1 G . (1 - 7i) V « rf - w*)) ^ds + err 5 , 



with 



|err 5 | < G|| (l Ge V « d - <)) g|U 1(Ge) ||VW £ - VWJ L - (G .) 
< C|G e |||HI^(n)l|VW e - VW*|| L » (Gl) . 
Note that thanks to Proposition |6j we know that 

||VW s -VW*|U~( Ge )= (l)p6|| G1 (r5), 
therefore |err 5 | = o(l) |G £ | P&|| G i(x5)|M|vK 2 .4(fi)- Finally, by Remark[7j 
1 f d f 

w~\ J vw * ■ (i ~ 7i) v ~ ^ ~^ dx = y vw * ■ (m * )T Vw dM + erre 

where 

|err 6 | = o{l)\\S7 W*\\ c o, 1/20 {a) \\S7w\\ c0 , 1,20 {S2e) = o(l) \\$b\\cHn)\\M\w 2 ^(n)- 



Combining (|62| , |63j) , (|64) , |65| , (|66j) and |67j), we obtain (|61J). □ 
Proof of Theorem [^} A straight-forward integration by parts shows that 

^s($b)= [ (71 - 1) VW e , d ■ VW e dx, 

which we rewrite in the form 

& e {$b)= ( (71 - 1) VW» • VW» dx 

+ I (a e , d - a e ) V (W Sld - W £ ) ■ VW* dx + r £ , 
with r e = (71 - 1) (VW e • \7W e - VVF* • \7W*)dx 

+ I (a e , d - a e ) V (W M - VK £ ) • V(W E - W,) dx. 



By Proposition |6j Lemma [2] Theorem [3] and ( 48 1 , we have 

\r e \ < C \G £ \ ||VTU £ - VW*||oo (ne) . \\$ b \\ci[ti) = o(\G e \). 
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On the other hand, Proposition [7] shows that 

{a e>d - a £ ) V {W £ , d - W £ ) • VW* dx = - \G e \ f M*VW* ■ VVF* dfi + o(|G e |) \\W4 w2 , Hn) 

and this establishes the representation formula given by Theorem |4j thanks to Theorem [3] with 

M ij = ( 7l -l)S ij -M* j . 



The bounds (52 1 on M* imply the announced bounds on M. □ 

To conclude this section, we now provide an alternative characterisation of M* , following |15j . 

Proposition 8 Let M* be the polarisation tensor introduced by Proposition^ Let ip be a uniformly 
positive, smooth function on Q , and ^ £ R 3 . Then M* satisfies 



(71 " 1} \e^, 



n 7i 



r—r mm 

|& e | ToeHj(O) 



a £ .d 



T1 — 1 

v™ + - IgA 

7i 



ipdx + o(l), 



where o(l) may depend on ip but goes to zero with e. 
Proof Let £ e be the solution in Hq (Q) of 

div (a M VVCe) = div (V (1 - 7i ) la J) . 
Note that Q is the unique minimizer of 



(69) 



Ti — 1 

V»+- 1g £ ? 

71 



•0 dx. 



In particular, 
and therefore, as 71 S C°(17), 



71 ~ 1 

7i 



(71 ~ 1} \e^ x 



7i 



(71-1) lel 2 



|£|^ + (1). 



/fi 71 

Let us prove that £ £ satisfies an estimate similar to that of Lemma [2] namely 

||vai iW + \G e \- 1/3 UeWmGe) < c \g e \ 1/2 , 



(70) 



(71) 
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where C may depend on ip, but is independent of e. Testing (69) against £ e , and integrating by 
parts, we obtain 

HVCe|| 2 LWi < J^J^a E ^j\V( E \ 2 dx 

1 ; (i- 7 i)e-vc £ 



mini/) J G 

<C\G E \ 1/2 \\VCe\\ LHn)3 , 

using the Cauchy-Schwarz inequality. On the other hand, using Holder's inequality followed by 
Poincare-Sobolev's inequality, 



\\Ce\\mG E ) < |G £ | 1/3 KeWmm <C|G £ | 1/J ||VCe||L 2(fi ), 

and (71 1 is established. Introducing V E = Y^=i( v l d ~ v l)£i' the identity (60l shows that 



,1/3 



a e4 \VV e \^dx= / {(l- 7l )l Ge )Zip-VV e dx + o(\G e 
n Jn 



= / O^VCe ■ VV E dx + o{\G E 

Jn 

On the other hand, using Proposition |6j followed by (71), used twice, 
a e , d \VQ 2 Tpdx= [ ((l- 7l )l G J^- VCecb 



We have obtained that 



= ^((l- 7l )l G J (X>"&1 ■VQiPdx + o(\G e \) 

= ((1 - 7 i)1gJ fe Vw &l •V(C^)«te + o(|G e 

= / a E ,^VK-VC £ ^ + o(|G £ |). 
Jn 

I a e , d \VV e \ 2 Hx= [ a e . d \VQ 2 ipdx + o{\G e \). 
Jn Jn 



The conclusion then follows directly from the above identity, ( 59 ) and ( 70 ) . □ 



5 Proof of the homogenization result 

This section provides a proof of Theorem [3j The proof is divided in three steps. 



First step: Derivation of the homogenization problem |13j ). 

This step uses the same ingredients of [12] and [5], but the boundary conditions are different. 
Choosing W E — Hq({2) as a test function in equation ([7]), we obtain 



/ a e \VW £ \ 2 dx = / a e VW E -V$ b dx 
' n Jn 



a E VW E ■ V<P b dx 



VW E ■ V<£ b dx. 
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This combined with (|3|), ( fl2"| ) and the Cauchy-Schwarz inequality yields 
/ a e \VW E \ 2 dx < c\\V$ b \\ L oo {Q) ( f a £ \VW £ \ 2 dx 



1/2 



hence the energy estimate 



/ a £ \VW £ \ 2 dx < c\\V<P b \\ 2 La 



■(«)■ 



(72) 



Estimate (72) implies that VW £ is bounded in L 2 {Q). Due to the Dirichlet boundary condition and 
the regularity of J?, the sequence W e is bounded in H 1 (S1), and thus converges weakly to some 
function W» in iJ 1 (J?) up to a subsequence. 

Integrating along vertical lines, we can prove (see [5], [12] for further details) that the rescaled 
function V £ = (nr 2 ) -1 1q c W £ converges weakly-* in VW(J?o) (i-e., in the sense of measures on i? ) to 
some function V* € H 1 ((— L, L); L 2 (wo)) • Moreover, the uniform repartition of the highly conductive 
cylinders Qm,n,e and the continuity of <Pb imply that V* inherits of the Dirichlet boundary condition 
on cj x {±L}. 



ndzV* weakly-* in AA(£2q), 



On the one hand, by ( 12 ) and <[72J) we have 
a £ d 3 W £ l Qc = d 3 (irr 2 a £ V £ ) 



and there is no transverse diffusion induced by the cylinders Q m ,n,e (see [H]). Moreover, due to the 
energy estimate (72) and the periodicity of a e , there is no concentration effect of a £ VW £ on dfl . 
Therefore, we obtain the convergences of the flux 

j a £ VW £ c% ->> VW» weakly-* in M{Q) 

| a e VW E (1 - c%) k9 3 K e 3 weakly-* in A^(/?o)- 



(73) 



On the other hand, the sequence c £ defined by (11) satisfies the convergences (see 



1 weakly in H 1 ^) and |Vc° 



7 weakly-* in At(i?o)- 



(74) 



Then, thanks to [HI Lemma 2] combined with ([3]), and again using the fact that there is no concen- 
tration effect on dfl , we get the convergence 



7 (W* - K) weakly-* in M(/2 Q ), 



(75) 



which induces the non-local effect in the homogenization process. 

Let i/jq, i>\ £ C\{fi). Choosing ip c° + i/>i (1 — cf) as a test function in equation (13), and passing 
to the limit as e — > thanks to the convergences (73), (74), (75), we obtain the equality 

VW, ■ Vip dx+ I nd 3 V* d 3 fa dx+ I 7 (W« - K) (V-o - fa) dx = 0, 



which corresponds to the weak formulation of problem (13). 
Second step: Proof of the local regularity of W* and V» . 



Rewrite the first and the second equations of problem ( 13 ) as 

- AW, = -7 loo (W* - V.) in O, 



and 



k d 2 3 V* + 7 K = 7 W* in J? = w x L). 



(76) 
(77) 
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As W* G H 1 {^2) and <P b € C l (Q), V* G i/ 1 (i7 )- By the Sobolev embedding theorem, the right- 
hand side of (76) is in L 6 (f2). Thus, as dfl is Lipschitz thus satisfies the exterior cone condition, 
the de Giorgi - Nash - Moser estimate (see e.g. Theorems 8.22, 8.27]), W* G C°' Q (J?) for some 
a G (0, 1). Consequently, V*, d X3 V*, d% 3 V* G C ' a (f2 ). 

Going back to (j76j), the right-hand side belongs to L°°(Q). Thus, W* G (f2) for any p > 2. 
Bootstrapping between (761) and ((77) we obtain G Cg c (f2 ) H Cj^ c (i7 \ J? ) and K G C^ C (V2 )- 



Third step: Proof of the corrector result (14). 

Denote by E s the left-hand side of (14). By equation |7]) we have 



a F \S7WA 2 dx = / a f VW f • V<2>h dx. 



Hence, by the definition (11) of cf, it follows that 



E e = / a s VW e -V$ b dx- 



|V<| 2 (W*-l4) 2 cfa; 



+ / «)^|VW*| 2 cfa; + / a e (l - c°Y (d 3 V*Y dx 



VW e ■ Vc a e {W* -V*)dx- 2 c% VW e • VW* dx 



!'->„ 



n 



a e {l-c%) d 3 W £ d 3 V*dx + o(l) 



Then, passing to the limit as s — > thanks to the convergences (73), (74), (75) combined with the 
continuity of the functions W* — V* and d 3 V*, we obtain 

lim sup E e < E a = l VW* • V<Z> 6 dx + K d 3 V* d 3 <P b dx 

e->0 Jn Jn n 



+ f(W,-V,ydx+ / \VW.\ 2 dx+ K (d 3 V*ydx 



O 

VW* ■ V<2> 6 dx 



j (W* - V*) 2 dx - 2 \VW^dx-2 n(d 3 V*ydx 



<? 



nd 3 V* d 3 <P b dx 



\VWJ 2 dx- 



7(W* - V*) 2 dx 



k (9 3 K) 2 dx. 



Here we have used 
a £ (l-c-) 2 - 



lim 



a e (l — c°) dx — lim (a e nr £ \ — k weakly-* in M(f} ). 



e-s-0 



On the other hand, choosing W 7 * — <$> b as test function in the first equation of ( 13 ) and V* — <l> b in 
the second equation of (13), it follows that 



VW, ■ V(W* - <P b ) dx + / 7 (W* - V*) (W* - # 6 ) dx = 
j.' Vr? 

« d 3 V, d 3 {V*-$ b )dx- j 7 (W» - V, ) (V, - <P h ) da; = 0. 



Therefore, adding the two previous equalities we obtain that E = 0, which gives the thesis. 
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6 On supremum estimates 

A less ad hoc version of Lemma [I] is given below. It is probably known to the experts. 

Lemma 3 For fl C W 1 , assume that w £ H^Q) satisfies 

-di(A ij \x)djw) + c(x) w<f + div(h) in fl , 

for some f € L p (f2), h <E L 2p (f2) n , p,p > ^. // the coefficients A 1 ^ and c are locally bounded and 
satisfy 



A 1 -* (x) £i £j > v |£| 2 for any x £ fi and £ £ 
c(x) > A for any x £ J? , 

2p— n 



then w + £ L°°(f2) and, for any q > 2 n-n an< ^ 9 ^ ' 



l|w + llL~(f2) < C{n,q,v) \ \\f\\ L pm) + -\- \\h\\ L 2 P{n) 
L A « A 2 « 

Proof of Lemma^ We follow de Giorgi's method for the proof. For the sake of Lemma [3] we will 
assume that 



/ a E f p dx < oo and / a E h 2p dx < oo 



for some p,p > 1; note that n = 2 in the present proof. 

For k > 0, let v4(fc) = {a; £ cj : ip e (x) > k}. Let £? represent the quantity 



Z? = 



A(fc) 



a e | V (</j e — fc) + | +A / a £ (<p s — k) + + A / a e fc (ip E — fc) + . 



A(fe) 



Integrating (22) by parts against (u — fc)+, we obtain 



lA(k) 

We shall use the notation conventions that 



A(k) 



B = I a e f(u-k)+- I a E h ■ V(u - k) + . 

A(k) 



(78) 



11/11- 



CP(A k ) 



l/p 



a s f p dx\ and N (A k ) = / a 6 da;. 



Using the weighted Sobolev embedding given by Lemma [4] we obtain 

Pe \\(<Pe - k) + \\ 2 C3(Ak) + A \\(Ve - fc)+||^ (Afc) 



< B, 



where 



Using Young's inequality, we have, for any 1 > 9\ > and 0% = 1 — #i, 



(79) 



> 



+ \\C 2 (A k ) 



26» 2 



+ ll£»(A fc ) l'^ £ fc )+H£ 2 (A fe ) 



(^-fc)f 



2/r 



(^-fc)f 2 



2/r 

C 2 'r°2(A k ) 
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On the other hand, Holder's inequality shows that 



+ \\C-(A k ) 



(<Pe ~ fc)+ r {fe - k) 



9 2 r 



2/r 

CHA k ) 



< 



(ip e - k) 



Sir 



for any 1 < a, (3 < oo such that ^ + 4 = 1. Choosing 



2/r 

£=(A fc ) 

a 



(fe - k) A 



2/r 



rO, 



2 

r0 2 



yields (3 = -§- provided 



1 _ 6»i 6> 2 
r ~ T + 2~' 



We have obtained 



^rii^-fc) 2 



(80) 



< || V (ip e k) + \\ 2 c2(Ak) + X \\{<p e k) + \\ 2 c2(Ak) < B. 



Let us now turn to the right-hand side. 
Case 1: h = 0. We have 

B< \\f(<p e -k) 



+ »C 1 (A k ) 



< \\f\\c P (.A k )\\(^ - k) + \\ cr{Ak) \\l\\ CHAk) 



with ^ + ^ + ^ = 1- We require that 61,62 be chosen so that k > 1. For p = 2 and s > 2, this 
requirement is fulfilled by any 6\ > since r <E (2, s). For the general case when p > §, we can satisfy 
the above requirement by selecting s sufficiently large (but smaller than the Sobolev exponent) so 
that - + - < 1. We have obtained 

p s 



life - k)_ 



< 



>+\\Cr(A k ) 

where x = r (l — P 1 — r -1 ) = r/z-c. Now, for all h < k, 

(w - h)~ 



f \\fh Hu) N(A k ) 



x/r 



(81) 



Therefore 



N(A k ) 



- fc )+| 



A k (k~ h) r 



dx = \\(ip £ - h)_ 



^ r ( A h) \k-h\ T 



C-(A k ) 



< 



1 



To proceed, we select 6\ such that 6\ > 
We now set kj =2d(l — 2~- 7 ~ 1 ). Introducing 



A J \k-h\x 
R^2)-Thcn X >l 



II/IUpmIK^-^+II^^) 



P = \\f\\&w(£) (y) 22Xd_x ' ^^ = ^11(^-^+11^) 



Then (81) takes the form 



x j+ i < y*x x „. 
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An induction shows that 

if x < 2- x/( - x ~ 1)2 7 then Xj < T 
Therefore, lim^oo Xj — 0, and — 2d. Thus, d is given by the constraint 

x < 2~(A % . 
Using the L°° bound just derived, for j = 0, we obtain 

Therefore we choose d to be 



Altogether, we have obtained 

02 



( Q \ /i 

^ c '(M>«+.X > » , )ll/ll J cp( tlJ ) ( y) ( 



For p = 2, we write 0j = ^ 2 {s-2) > with 1 < /i < 2i£_21 < 2 and use (79 1 to obtain 



||^|| L oc. M < C(fl, S) \\f\\ C 2 (uj) — gg— . 

Now, for a fixed 1 < a < 2, we can choose s > 2 sufficient large so that we can choose /x = a 
in the above estimate. By further enlarging s, we see that we have obtained the assertion for /3 as 
close to 1 — % as we wish. The conclusion for smaller j3 also follows. 
Case 2: / = 0. Using Holder's inequality we get 

b < Uhiu^ ||v(u - fc)+IUw)) ^ (Mk)) 1 ^ ■ 



Together with ( 80 1 , this shows that 

0i/2 / „ x e 2 /2 



\\(^- k ) + \\c W ^ (^) 1 (x) 2 l|h||^ M ^V(A(fc)) x/r , (82) 

where ^ = | — Note that (82) is a reverse Holder inequality of the same type as (81). Therefore, 
arguing as above, and provided g — ^ > £ , to ensure that x > lj we obtain 



||V»e|lx,»(o;) ^ C(\u\,K+,X,r) \\H^P( U ) (j^J ( 



0l/2 / n \ e 2 /2 
2 



A 

Considering now the case p = oo, introducing 6\ = (a — 1) 2 / s s _ 2 \ , with 1 < a < 2^^, and /3 — 
we obtain 

II^IL-H < ll h IL«( W ) £ (a-l)/2 A ff ' 

as announced. □ 

Proof of Lemma^ The proof is similar to the one above, with a E replaced by v. The constant fi(s) 
is now independent of e. Some of the details are provided in the proof of Lemma [T] for the reader's 
convenience. □ 
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Lemma 4 Assume that conditions ([I]), ([2| and ([3| hold. Then, for any s € [2, oo), i/iere exists a 
constant C(s) > smc/i that 



a £ \v\ s dx <C(s)e 2 - s [ / a £ |Vw| 2 dx V« € ^(w). 



(83) 



Proof As a first step of the proof, we establish a rescaled version of (83 1 in the cell Y = (—1/2, 1/2) 2 . 
Define a rescaled conductivity A e (y) = a e (ey), for y G Y. Using a r e -rescaling and the Sobolev 
embedding of H 1 (D 2 ) into L S (D 2 ) applied to the unit disk D 2 , we get that 



yveH\Y), 



JD{r e ) JL 



V-+ V 

D(r e ) 



(I 



dy<Cr s E [f \\7V\ 2 dy 

>D(r e ) 



Moreover, by estimate (3.13) in jTTj we have 
We H\Y), 



V - 



V 



D(r e ) 



dy < C|lnr e |5 / \VV\ 2 dy 



Then, combining the two previous estimates it follows that for any V € H 1 ^), 



1 


V -Jr V 


s 

dy<C 


[ Vdy-S Vdy 


w 


V-f V 


)D(r e ) 






■Jy JD(r E ) 




JD(r € ) 



dy 



<<7|lnr e |* ( J \VV\Uy) +Cr s e 



j 



\VV\ 2 dy 



Hence, by the definition of A £ and ([!]), 

' dy<C\lnr £ \i ^A £ \VV\ 2 dy 
Since the Sobolev embedding inequality in Y gives 





v ~L v 


lD(r E ) 





I 


V -L V 


'*<o(l 









dy<C[ / \VV\ 2 dy <C( A £ \VV\ 2 dy , 



we thus deduce the following estimate in Y, 



VVeH l (Y), / A £ 



I- 


v ~L v 







dy<C\\nr e \-* / A £ \VV\ 2 dy 



(84) 



We now turn to the proof of ( p3| . Let v G Hq(lj) and extend it by zero outside w. Rescaling 
estimate (84) in each square e(n + Y), for n £ Z 2 , we obtain that 



a £ v - v £ 



'dx < Ce 2 |lnr e |3 ^ ( / a £ \Vv\ 2 dx 



e(n+Y) 



< C e A \1nr s \* [ / a £ \Vv\ 2 dx) (since s > 2), 
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where v e is the piecewise function which takes the average value of v in each square e(n + Y). 
Moreover, we have 



< 

Tie 
< C 



I / a £ dx ) 


+ vdx 


\Je(n+Y) J 


Je(n+Y) 



._™ 9 Jei 



\v\ s dx (by the Jensen inequality) 



s{n+Y) 

\v\ s dx 

< C [J \Vv\ 2 dxj 2 (by the Sobolev embedding of H%(u>) into L s (uj)). 
Finally, combining the two previous inequalities we find 



'dx < C(e 2 I In; 



1 



a F Wv\ dx 



which yields the desired estimate ( 83 ) taking into account ([I]) . □ 

Remark 8 The constant in the weighted inequality of |33j provides an estimate from above of the 
constant appearing in (83). It is not clear that this constant is optimal. The dependence in e of the 
constant in (83) is optimal: this can be verified with the choice v = g e , where g e is given by (27). 



7 Proof of Proposition [2] 

All of this section is in the two-dimensional setting. We will therefore drop the subscript 2 to denote 
two-dimensional gradients or divergences. We consider the solution u e of (21), for A > Ao > 0, 

J — div(a £ Vti £ ) + X a £ u £ = f + a e g + div(h) in to , 
1 u e = on dui . 

Let us start with a simple energy bound. 
Lemma 5 We have 



X\u e 



dx < 



[l/l 2 



\g\ 2 ]dx + C / \h\ 2 dx 



Proof Integrating (21) by parts against u e yields 

|Vu £ | 2 + A|w £ | 2 dx 



fu £ + a £ gu e — h • Vit e 



dx 



(85) 



Using the Cauchy-Schwarz inequality on the right-hand side we obtain ( 85 1 . □ 

We now turn to the main part of our estimate. Our strategy, inspired by the limit case when 
the conductivity tends to infinity independently of the periodic structure [18,8,16!, is to consider 
three contributions to u e . The first one is the contribution of the right-hand side when no highly 
conducting fibres are present. We introduce v £ Hq(lj), the solution to 



Av + Xv = f + g + div(h). 



(86) 
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The second is the contribution coming from the average of u £ — v on the cross section of the fibres. 
Set 

Um - n = Tan 1 / [u e (x) - v(x)]da(x), 

\0-L> m ,n,s\ JdD m ^, e 

and define u £ £ \ D s ) as the solution to 

—Au £ + A u £ — inw\ D £ , 

u e = on du , (87) 

&e ^m,n On dDrn,n,e • 

The third is simply the remainder, given by u £ — u £ — v — u £ on w \ D £ . We will consider each 
contribution separately. Intuitively, our estimates show that, outside the fibres, u £ consists "mainly" 
of two parts: v, the contribution of the background medium, and u £ , the contribution from the high 
contrast. 

Proof of Proposition^ Fix r > K£ 2 < 1+r '> . 

Since u £ + u £ = u £ — v, Lemma[7] (with (3 — n/2) and Corollary [4] show that 

\\u £ - w|| L ~( w =) + ||Vu e - Vu|| L oo (<) 
< ( y (Il/IU 2 H + |lv^slU 2 M + ll h IU~M) + IK - v\ L <*> {uj >) 1 

where cjo <e <e Using the L°° estimates on w given by Lemma [6] (with q — 2/77) , we obtain 
IK - v |U°° (w? ) + ||Vu £ - Vu|| L oo (<) 
C 



< 



I (l^H i2 ( w ' + HV^SlU^o;) + ||h|| iD o (a)) ) + llttellioa^/jj 



Lemma 10 together with Corollary [4] show that provided v < 77, 



C e 2 _ 

[V 2 u e - V2u]c(w£) < tuttIK ~ w|U«. (w /) + C-3- (||/||l2 (w) + HvWfflU^u) + ||h||ioc( w) ) . 

A 2 2 ,\ 2 

Using again the L°° estimates on v given by Lemma [6j we obtain 

[V 2 u e - V 2 u]oKQ < ^rrv- f^a (ll/IU 2 (u>) + II x/aT^lli^^) + ||h|| £ «. (w) ) + |K|| i0 o ((</) ^ , 
which concludes the proof. □ 

7.1 Estimates for v 

It is quite straightforward to obtain estimates on v, as the following Lemma shows. 
Lemma 6 There holds 

V / A||V-y|| L 2 (tj) + A \\v\\ L 2 (uj) < C \\f\\ L '(ui) + IMIz^w) + V^ll h IU 2 (u;) , 
< C{q) [||/||i2 (w ) + ||3||l2( w ) + ||h|| LO o (aj )] , 2 < q < 00. 
Furthermore, for any q > 2, and any uj 1 <s uj, 

IMU- ( -o < [II/IIl2 M + ||<?II^m + IWU-m] ■ 

A? 
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Proof The L 2 estimate for Vi> and v follows directly by testing (86) against v. To obtain the L q 
gradient estimate, we write v = v\ + V2 where Ui,«2 € Hq(u:) are the solutions of 

— Av\ = Xv + f + g = f and — Av2 = div(h). 
Next, by standard elliptic estimates, we have 

WviWm^) < C||/||i2 (w ) < C X\\v\\ L 2 (u]) + \\f + gWwfa) < C (||/||i,2( w ) + ||flr|U a (w)) , 

IIHIw^M < C INU«(w)- 

Since H 2 (lj) «-> W 1,q (uj), the second estimate in the lemma follows. The last assertion follows from 
Lemma [H □ 

A direct consequence of the above result is a local L 2 bound on Vv. 
Corollary 3 For any u' <g cj, f3 < h and g < 1 ~?^ > , there holds 



\\Vv\\ LHul/) <C\u>'\ 



(ll/IU^) + \\V^9\\lHu)) + l|h|U~M 



Proof We split v — v a + Vb where V a ,V}, € Hq (uj) are the solutions to 

— Av a + Xv a = f + g and — Av\, + Xvb = div(h). 
By the L q gradient bound in Lemma [6] and the Sobolev embedding theorem, 

I! VV a \\LP <C{\\f\\ L * ) and || VvbWLP^ < C \\h\\ L ^(uj) for any p > 2. 

Thus, by Holder's inequality, we have, for any ui' C uj and S £ (0, 1), 

||Vv„||i» (u0 < C (II/IIl^m + Ibll^Cu,)) Kl 1 " 5 and || V^IU^j < C||h|| L » M la/1 1 - 5 . 
On the other hand, by Lemma [6j we also have 

Q 

||Vu ||i3( w ) < ^ (H/IU*(uO + IMU^)) and || Vu 6 || L 2 M < C ||h|| L2(w) . 
The conclusion follows by an interpolation. □ 



7.2 Estimates for u s 

The first part of our estimate for u £ is given by the following Lemma. The second part, concerning 
the Holder regularity of its gradient, is given by Lemma [l0| 



Lemma 7 For any < (3 < k and r > k£ 1 + 2 ' 3 , the solution u e of (87) satisfies 



C 



l|We|U°°«) + ^ 1/2 ||VWe||i~( w? ) < SUp |u m ,„| , 



where the constant C depends on (3, k and to only, and uj\ is defined in 

Here and below sup m n refers to the supremum taken as (m, n) varies in I £ . 

We shall use two local estimations. We first estimate the gradient Vu e on the outer boundary du>. 
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Lemma 8 There holds 



Proof By the maximum principle, 



< C Sup \u min \ 



SUP \U S \ < SUp|M mi „| . 

uj\D e m,n 

Fix wo ^ w' g w" g w. Let ?/ be a cut-off function which is one in uj \ u>' and vanishes in ujq. 
Using r] 2 u e as a test function in ( 87 ) , we obtain 



?7 2 |Vu e | 2 + 



Xrful = -2 



Thanks to 



we conclude that 



rju^rj ■ Vu e 



< 



7 ? 2 |Vu £ | 2 + 2 



u 2 |V„| S 



l|Vu e || L 2 (w \ w , ) < C sup|u mi „| 

m,n 

Let € \ u/) be the solution to 

-zkt* =0 in u \ uj' , 



(89) 



,± - 



u = 



on dui , 
on eta/ 



where uf — max(±u e ,0). By i? 2 -estimates, valid up the boundary dui thanks to the vanishing 
boundary condition and the regularity of duj, u ± £ H 2 (uj \ oj"), with 

\\ u± \\h 2 (uj\uj") < ^H&ellfl-iOAa,'). 

This, in turn, because satisfies a homogeneous equation in uj \ uj' , implies a better estimate, 
namely a VF 2,p -estimates for an arbitrary p [2Q., Lemma 9.16]. For p large enough, thanks to the 
Sobolev embedding Theorem, u± € W 2 ' p (uj \ uj") ^ W 1 ' 00 ^ \ uj") and we deduce that 



"\<C\ 



u e \\m(uj\u>>) 



Inserting (88) and (89) in this last estimate, we obtain 

Hit* || wi.°°0A<"") - Csup|M mj „| . 

As an application of the maximum principle we note that u > and so 

-Au + + Xu + > = -Au s + Xu £ > -A(-u~) + X(-u~) in u \ uj' 
Applying the maximum principle again, we thus get 

u + > u £ > —u~ in uj\u)' . 



(90) 



Since the three functions agree on duj, the required estimate follows from (90). □ 



Next, we estimate the trace of the gradient Wu e on the boundary of each rod D n 
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Lemma 9 There holds 



where 



|Vw E ||Lo°(aD m , n , s ) < C6(X,e) swp\u r , 



0(\,e) = < 



£ 



— -, ^ if V A = with < a < 1, 

r E 1 - a + e 2 



1 



(91) 



— < V A. 

er F 



This result is proved in Appendix |A. 3 



We are now in position to prove Lemma [7| Note that the fundamental solution <P of — A + A 
on K 2 is given by 



$(x)=*(\x\) = ±K (V\\x\), 



(92) 



where Kq is the second modified Bessel function of order zero. In what follows we will make extensive 
use of well known facts concerning K and its derivatives, easily recovered from the recurrence 
relations and asymptotic properties, see e.g. (32] ■ Let Kq denote de n-th derivative of Kq, n > 0. 
Then |1£q | is decreasing, and, for all x > and a > 0, 



\K (x)\ < C(a)^ and |**0(*)| < -^^,n 1 I 



1 + x° 



x n (l + x a )'' 



(93) 



Proof of LemmaVA Our argument uses layer potentials. We set M = sup m n \u m<n \. Let # be given 
by ([92]). We have 



u e (x 



d(u\D e ) 



]^f^(y)®( x ~y)~u £ (y)J^(x-y) da{y) 



(94) 



Here n(y) denotes the outward unit normal to d(w \ D E ) at y. Fix any multi-index J £ Z 2 . Differ- 
entiating (94 1 we obtain 



W J u £ (x) = 



d(w\D E ) 
T J _|_ T J _1_ J J 



^(y) viM? -y)- uM ^(x- y )\ My) 



(95) 



where 



du £ 
dn{y) 

(m,n)e/ e 



(y) V J x ${x-y) da(y), 



d\7 J $ 

, ll! n My) d^) iX - y)d(j{y) ' 



dD„ 



dn(y) 



(y) V J x ${x - y) da{y) 
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Fix x 6 <jj% . We now estimate the right-hand side of ( 95 1 . We begin with the easy estimates for 
l( and 1%. Thanks to (931 and Lemma pM 



M for all $ > 0. 



(96) 



To bound Jj 7 , we estimate the contributions of inclusions located "radially" between er + je and 
st + (j + l)e for j = 0, 1, . . . We obtain using (93), 



Ce~ 



\H\ < C J2 M (j + 1) K^ +1 \VX (j + r)e) 

Ce- 1 



er F 



< CM 

< CM 
= CM 



\P T \j\+i+2p £ \j\+2p 



71+2I E (^Cj + r)e: 



|J| + l+2/3 



3=0 



-. ^ , sup t\ J \ +1+ ^K [ ^ +1 \t) 

X T \J\ + 1+2P £ \,J\+l+2f} 0<t ^ OQ W 

/? > 0. 



(97) 



^/3 T |J| + l+2/3 £ |J| + l+2/3' 

Let us now turn to the third integral, Using Lemma [9] and counting contributions from "rings 
of inclusions as above, we have 



\li\<Cj2 0(X,e)M(j + l)V\ lJl \4 lJl \VX(j + r)s) 
3=0 



er R . 



Since #(A, e) can be "large" , we proceed to decompose // into two parts: one counts the contribution 
from the ring of inclusions closest to x, where the dependence on t is dominant, and the other counts 
the contribution of the further away inclusions, where r does not play a role. Proceeding in this way, 
we compute, using the decay properties of Kq for the second term, 



|#l < C ]T 0(X, e) M (j + 1) V\ lJl ^ |J|) (\/A (j + r)e) 



3=0 



<CM9{X,e)er e Vx lJl K { lJ{) {VXt e) 



X J — near neighbour terms, i.e. j = 



CM8(X,e)er e jVx K ( lJl \VXje 



Y J — distant neighbour terms 



The near term X can be bounded as follows, 



X J < CM 

< CM 

< CM 



- 9{X 'f )r£ .(VXre)^ 4 lJl \VXre) 

\0 T \J\ +2/3 £ \J\- 1+2/3 y > \ J 

0(X,e)r £ 



\/3 T \J\+2/3 £ \J\- 1+2/3 t> % 

0(X,e)r £ 

\0 T \J\+2P £ \J\-1+2P' 



\J\+2P 



4 lJl) (t) 



p > 0. 



(98) 



:i8 



Similarly, we have 



Y J <CM g(A ' £)rE . Y (VXje)^ K^Wlje) 

~ \0 r J -1+20 



A /3 £ |J|- 1+2/3 



<CM g(A ' £)r l suptl J l+ 2 ^ 

A /3 £ |./|+2/3 f> £ 

<CM^, /3>0. 

A /3 £ |J|+20' 



(99) 



We then insert estimate (91) in (98) and (99). For A > — , we can choose (3 = 1 and obtain 



< CM 



< CM 



For A < — we can take j3 = \ to get 



(100) 



(101) 



([lOOj) and |101J ), we obtain, for all A > A 

CM 



These estimates are enough to conclude for the case J = (0, 0). Adding (96), ( |97| ) (with fi = |), 

A > A , 

. ^CMf_. r F e r r 1/2 / * 
tie (a;) I < 



X A 



1 



which proves the first desired estimate as r > ke 1 + 2 ' 3 ^> e. 

For | J| = 1 and A < estimate (101) deteriorates. We keep the near estimate (98) for X J and 
improve the far estimate (99) for Y J . First, we write 

Ce- 1 

y (o,i) + y(i,o) =CM6(\,s)sr e ^ jVx K^\VXje) 

3=1 



< CAf 



Note that since t — > t\K^p(t)\ is decreasing and summable on (0,oo), its lower Riemann sum is 
bounded from above by the continuous integral. Thus 

Y (o,i) + y(i,o) < c M 6{\e)r e t°° ^ R ( X) ( ^ ^ 

6 Jo 
9(X,e)r e 



CM 
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This inequality together with (91 1 and (98) yields that for A < — 



\i^\ + \i^\<™( e ^Z + i), o<p<l 

13 1 13 \P V T l+20 /' 2 



(102) 



Choosing (3 = {3 in the above estimate and in (96), (97) and recalling (100), we conclude for all 
A > A that 

, \ i CM ( r £ e 1 - 2 ? rl /2 / e \\ 

\Vu E (x)\ < —\i + T2+2[j ^ +2f} + _ w + — + 

1-2/3 

which provides the remaining estimate as r > Ke 1 + 2 f . □ 

We now turn to the derivation of Holder estimates for the gradient of u e . 

1 1-2,3 

Lemma 10 For any Q<j3<^,Q<v<2[3,n>Q and r > keW-i there holds 

C 

[Vtt e ]c«) < Tfl^iT sup|u m! „| 

A 2 m,n 

/or some constant C independent of e, r and A. 
This Lemma is proved in Appendix | A. 4 1 

7.3 Estimates for ii t 

The residual term u F satisfies 



—Au e + A u e — in cj \ Z3 e , 

u e — u £ - v - u m>n on dD m<n>e , (m, n) e I e 



(103) 



We derive the following estimate, which shows that u e is negligible compared to other contributions 
to u r . 



Lemma 11 For any J £ I 2 , r\ > 0, and t > ne \ J \+ 2 i r 6 16( J|+2,;) ; the solution u e of ( 103 ) satisfies 

Crj 



\Vu. 



e\\L°°(u% ) 



< 



l/IU 2 (w) 



where the constant C depends on J , T], k and U) only, and ijj e T is defined in Q . 

For the proof of Proposition [2j we shall use the following simple consequence of Lemma [TTj 



Corollary 4 For any r > ke , the solution u £ of (103) satisfies 



|We|U>»(wf) + ||Vu £ || L o 



(«?) + [V« £ Jc».i(^) S — 
A 



< 



Ce 2 



— [||./1U 2 M + \\V^9\\lHu) + ||h|| L oc (w) ] 



where C > depends on k and oj only, and w% is defined in ^ . 

2 a 

Remark 9 A similar statement is possible with replaced by with arbitrary a > and /? > 0. 

However, we chose those powers to fix ideas, since, for example, IMIl 00 ^) an d ||'Se||i«>(a) E ) have 
slower decay in A, and the restriction on r (with respect to e) is more stringent. 



The proof will use the following result, whose proof will be given later. 
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Lemma 12 Given a G H 1 / 2 (S 1 ), and s > 0, there exists a unique solution in H 1 ^ 2 \ 5(0, 1)) on 
the exterior of the unit disk of R 2 of 

-Au + su = in R 2 \ 5(0, 1) u = aon§ 1 , 

still denoted by a, which decays at infinity. It satisfies 

\a(x)\ < Ce-^l/ 4 \\a\\ L 2 (sl) for \x\ > 3. (104) 

Remark 10 We do not claim that this estimate is in any way optimal, but is it sufficient for our 
purpose. 



Proof of Lemma 11 We first derive a bound for u £ away from D £ . For a given (m, n) £ 7 e , let a+ „ 
be the solution to 

"H» + Xa ti,n =0 m M 2 \ An,n,s , 

on dD n 



a m nC 2 -') ^0 as — > oo . 



By the maximum principle, a+ n is non- negative in to \ D £ . Thus, again by the maximum principle, 

(105) 



u £ < a m,« inuj\D £ 

(m,n)6/ e 



Now, thanks to Lemma 12 wc have 



|a+ | < C exp - 



| u e - v - u m<n \\ L 2( dDm n ^ in cj t/4 



On the other hand, by Poincare's inequality, (85) and Corollary |3j 

||u E -v- u m ^ n \\ L 2 {dDmn ^ < C y/er^ ||Vu e ||z,3(D miraiE ) + || Vu|| L 2 (Dm ^ e) 

< C^Vl (r £ + (er £ ) 2e ) x 



X" 



(ll/IU 2 (^) + 117% .glU^)) + l|h|U~M 



for (3 = j£ (0, an d £> = S (0, 1 4 )■ Inserting this estimate in the upper bound of a+ „ gives 



"L < C exp 



VAST 
16~~ 



(r £ + (er £ ) 2 e) 



Substituting the above estimate into (|105|) results in 
u £ < C exp 



^Xer\ {r £ + (sr £ ) 2 e) £ 



16 ; e 2 

A similar lower bound is derived by repeating the argument. We have obtained 

VXeA (r £ + (er £ ) 2 e) 



\ue\\L^{^ ,.) < C exp 



10 



(106) 
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We now proceed to prove the required derivative estimate for u e . 

We will only bound £ = d Xl ii e . The other (higher) partial derivatives can be bounded similarly. 
Fix x £ i>j%. Note that the disk B 3eT / 4 (x) C <^/ 4 - Using a suitable cut-off function which is one in 
B £T /2(x) and zero outside B 3eT /^{x) and a standard energy estimate, it is easy to show that 



C 

\^U e \\ L 2 {BcT/2{x)) < — ||u e || L -- V>v /;.. ,,. 



which together with ( 106 1 implies 



|Vu E || L 2( Be , 2(x )) < C exp 



V\er \ (r £ + {er e )^) 
16 " / e 2 ; ' 



(107) 



On the other hand, by differentiating (103), 

-Z\£ + A£ = in uj\D e . 
Usual interior De Giorgi estimates then show that 

C 

H£lU°°(B ST/4 0r)) < —U\\L2(B ET/2 (x))- 



Combining this bound with (107) and noting that x is arbitrary in o>£, we get 



ICIU-K) < C exp - 



VXer\ (r E + (Er e ) 



2o\ 



16 



The required estimate for £ then follows from the simple inequality e l x l < . Further higher 
derivative estimates are obtained by repeating the above process. □ 



To conclude, we now provide the proof of Lemma 12 



Proof of Lemma 12 It is straightforward to check that if a € L 2 (S 1 ) has the Fourier expansion 

oo 

a{9) = [afc cos fc# + b k sin fc#] , 



fc=0 



then it extends to a solution of (—A + s)a = on R 2 \ D\ which vanishes at infinity by 

oo 1 

a(r, 0) = T ^ - \a k cos kO + b k sin kd] K k (^/s r) , 

where K k is the second modified Bessel Function of order k. The function K k has a representation 
as follows (see e.g. [32]) 

/"OO 

iC fe (i) = / e - tcosh ^ cosh/crdr 



This implies the obvious bound 



A; > 0. 
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On the other hand, for k > 1, we also have 

K k (t) <£e-i*>* e- dr = * J~ e - dr = , 

and 

It follows that 

^fc(\/s) - r fc 

This is also valid for fc = 0. 

In particular, the above two bounds show that the series for a(r, 9) converges absolutely for r > 3 
and, 

00 ofc/2 

Hr,0)| < Ce-W-W ^J2 Q°*l + IM1 

fc=0 

< || a || ia(S1) forr>3 (108) 

as announced. □ 



8 Interior estimates for the boundary value problem 

In this section, we detail regularity estimates for the solutions W Ej i and W Ei 2 of the boundary value 
problems ( 30 1 which appear in the proof of Theorem [5] 



8.1 Regularity of W E ,i 
Our result is the following 

1 — ri 

Proposition 9 Let Q% = u% x (-L,L). For k > 0, r > ke 2 ( 1+ '» with 77 <E (|,l), and < v < 



2 (77 — |j, t/ie solution W E) i of (35) satisfie 

\\We,i - Vi\\ LO o m) + ll VT ^.i - V ^IL-(fl ? ) ^ Cr(«,»?)ll*fc||ox(/j), 
[VW' e ,i-V7i] c „ (n . ) <C(« J » ?> i/)||# 6 || C i (fl) . 

Proof We apply a Fourier series decomposition in S3. We write 



,X3) ~ 2^ ^-M,™^ ) Sm (,~2~Vy + 1 



Vi(a;',X3) - ^2vi, n (x') sin (^"(^ + l' 
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As W e> i — V\ is harmonic therefore smooth in 17 \ Q £ , the Fourier expansions of W e .i — V\ and 
V(W E) i — Vi) converge pointwise to themselves in Q%. 
Now w £} i, n € Hq(u>) satisfies 

2 9 

div 2 (a e V2lZ) e ,l,n) + ~^j2~ a eW e ,l,n = dfv 2 (h n ) 

where h„ = [h\, h%) is given by 

K = fwi^(~ (f + !))■ 

Since 933(^1 = 0, we see from an integration by parts that 

IKIIl~m < £||Wi|U~(fl) < ^WMcnoy (109) 
Lemma [3j applied with p = 00, shows that for any 77 < 1, 

C 

< _||j£|| ioo(u) , (110) 
Proposition [2] together with (109) and (110) now shows that, for every < v < 77 < 1, 



C 

||* E ,l,n - Ul,n|U°°(i<) + ||V2W e ,l,n ~ V2«1,7i||l°°(u;?) < 2 7;+l ll^llcH")' 

and 

C 

[V 2 u) e ,i,„ - V 2 u]c-K) < ri27)+1 _J l^fe|lcH^)' 

Following a variant of the proof of Theorem [I] we obtain consecutively the following estimates using 
the indicated choice of 77, 

- ViLoo^.) + ||V 2 W E ,i - V a Vi|| £a<>( „ f) < C||* 6 || ol(fl) with any 77 > 0, 



v 



[V 2 We,i - V 2 Vl]a<'(^) < Cll^bllc^r?) with an Y V > ^ ' 

H^We.i -^3^i|| i0 o (fi?) < Cll^bllc^fi) with an y ? 7 > 

[^3 1^,1 - d X3 Vi]c{n^) < C\\$ b \\cHfi) with an ^ r ? > 2 + \ m 

This completes the proof. □ 

8.2 Regularity of W e , 2 
Our result is the following 

1 — 7] 

Proposition 10 Let Q% = u)% x (-L.L). For n > 0, r > K£ 2 "+'i ttj^/i 77 € < v < 

3(77— |), i/ie solutions W e _2 andV e ^ of (38) and (43), respectively, satisfy 



\\w E , 2 - v s , 2 \\ LOO{m) + ||w £ , 2 - w £ , 2 |j LOO(r2?) < c( K ,77) e CT ||<? fc || cl( ^, 

[VVF £ , 2 - VK, 2 ] c „ (fif) < C(«, 77, v)e°\\$ b \\ 
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Proof As before, we decompose W £t2 in a Fourier series along the third direction, i.e. 

oo 

71 = 1 
71=1 

As usual, Wg, 2 — V e _2 is harmonic and regular in Q\Q e . Thus, we can sum the estimates on the 
Fourier coefficients to obtain estimates on W £ , 2 — V e ,2 and V(W e ,2 ~ Ve.i) in 
The problem satisfied by w e ^,n is now, in Hq(lo), 



-div 2 (a e V 2 w £ . 2 ,n) + \a £ w ea , n = a £ A £in + a £ B £ . n + C e , n , 



2 2 

where A = r -jfr, and 



A s , n = J ^sin (y\(t + L)) d 3 C(;t)d 3 W e , 2 (;t)dt, 
B s ,n = J sm(y\(t + L)^d 33 Ci;t)W e , 2 (;t)dt, 
C e , n = y sin + t)di. 



Proposition [2] together with Proposition [TT] now shows that for all < v < r\ < 1, 

C 

\\Ws,2,n ~ U2,n|U°°(w?) + ||V 2 W £ ,2,n - V 2 U2,n||z,«>(a>f) < £ CT -3^ 1 1 ^6 1 1 C 1 (fl) ) 

and 

C 

[V 2 W e , 2 ,„ ~ V 2 U2,r i ]c«) < e<T ^3^;l!^llci(r2)- 

Following a variant of the proof of Theorem [T] we obtain consecutively the following estimates using 
the indicated choice of 77, 

\\W e , 2 - ^2||i« (nf) + ||V 2 VF £ , 2 - V 2 V2|| ioo(fl?) < C^ll^bllcHfi) with an y f > 

[V 2 W £ , 2 - V 2 ^a]c"(n«) < Ce' 7 H^llc^n) with an y ^ > 3 + |> 

11 - ~ 11 2 

1)^3^,2 - 9a; s ^2 || ioo( „ f) < Ce^Pbllc^fi) with any 7? > -, 

[9x3^,2 - d X3 V 2 ]c»(n' T ) < C'e <T ||<Z'fc|| c i (f 5) with any 77 > ? + |, 

which concludes our argument. □ 
Proposition 11 We /iaue £/ie following estimate 

C 

||v^^e,n|U 2 (a;) + II V^e^n \\ L 2 (u) + Hv^^nll^^) < y^"^ ll^ll C?i (S2) ■ 
As a consequence, 

\\u> s ,2,n\\ L ^( U ) < C{ri) — \\$\\ C1{D) . 

for any i] < 1. 
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Proof Note that (42) shows that C £ . n satisfies 



C 

||C £ ,n|U°°(u;) < -j= \\<h\\L<*>(ii) 



Turning to B £jn an integration by parts shows that 

C 



B e ,n\ < I ||W e ,2|Uo. ( „) + / l U3X{ -l>,l>)\d 3 W ei2 (;t)\dt 



L 

-L 



which in turn shows that 



2 j„l - C 



a £ \B s ,n\ dx < 



\W Ei 2\\l<x>({j-) + 



a £ \d 3 W £ , 2 \ 2 dx 



UJ 3 X ( — (',(') 



Thanks to (39 1 and the interior estimate (41), we have obtained that 



C 



\\V^s B e,n\\L 2 (u) < ~y= II ^2 || L oo • 

We proceed with A £>n in a similar way. After an integration by parts, we derive that 



< 



C 



Whv) (|W e , 2 (-,0l 2 + \d 33 W £ , 2 (;t)\ 2 ) dt 



Multiplying this quantity by a £ , and integrating we obtain thanks to (39) and (41) 



Q 

a £ A £ . n \\ L 2 {LJ ) < H02|| L oo (f2 ) 



We have obtained that 



C 



W\faeA £ . n \\ L 2( U) ) + \\^fa £ B £ ^ n \\ L 2( u - ) + || C E , n ||z,2 ( w ) < — H^allioo^ • 

The first assertion now follows from Proposition [5] 

Finally, using Lemma[TJ we know that for any 1 < a < 2, and any j3 < 1 — a/2, 

\\We.2,n\\ L °°(Lo) ^ (llv^^,nlU 2 (o;) + II V®e B e,n II £=>(«) + ||<? e ,n|U=») 

<C(a,/3)-j^-||#|| cl(n) 
<C09) e *A- 1 / 2 -^||#|| OI(n) . 

□ 
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A Proofs of technical lemmas and propositions 

A.l Proof of Proposition [T] 
Note that we have 

J - AW* + 7 (W* —V*)=0 in Q , 
\ - k <9f 3 V; + 7 (V* - Wj) = in Q . 

And a simple induction shows that W* and V* are infinitely differentiable with respect to x% in J?o- 
To gain analyticity, we resort to estimates. Fix a ball B(p°,R) in Qq. We first derive an integral 
representation for W*(p°). For any r > 0, the function 

cosh(77|x-p°|) sinh^x - p°\) 

G p o (x) := — ^ Q| coth y^r) — ^ - 

4n\x — p u \ Att\x — p"\ 

satisfies 

(—A + j)G p o jT = S p o in B(p°,r), 
G p o )r = ' on dB(p°,r). 

A direct application of Green's formula then yields that for < r < R, there holds 

W.(p°)=7 / G p o, r (x)K(z)dz+- / W*(x)da(x). 

■/b( p V) 4?rr smh( v /7r) J dB ( p °,r) 



Here we have used 



dB( P °,r) 4irr sinh(y / 7r) 



It follows that 



W*(p ) / r smh(y/jr) dr 
Jo 

= 7 /' r sinh(- v /7r) dr / G p o !r (x) K(ie) ^ + 7— / W*(x)dx. 
Applying this identity to dx^ +1 ^W* we obtain 

,-R 

0& +1) W.(p°) / rsmh(^r)dr =7/ r sinh(^r) dr / G p(V (x) d(" +1) V; (x) dx 

JO Jo Jb( P °,t) 
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B(p°,R) 



7/ r sinh(v^yr) dr / G p o. r (x) d(" +1) V; (x) dx 

JB(p°,r) 

47r JdB(pO,R) \X-P0\ 
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Where we used the notation p° — (p" , p§ , ) . Noting that < G p o r (x) < cosh ^^2|^| p - , we deduce 



that 



d ( £ +1) W«{p°) / r sinh(77r) dr 



< C 



r smh(^r)dr / s coshes) ds\\d£ +1) V* \\l°°(b( p °,r)) 
Jo 



and we have bounded di™ +1 ^W*(p°) in terms of lower derivatives by 

C +1) W,(/)| < c[B 2 \\dil +1 W4L^(B( p o.R)) + ^ll^^ll^^o.S))]. (HI) 

Let us now turn to 14 . For a; 3 6 {p% — R, p® + R) we have 



5i: ) K(p°,p°,x 3 )=a + exp 



X3 



a_ exp 



%3 



where the constants a± satisfy 



\a±\ <C 



R 



\\d^v4 L ^ {Bip o. R)) + R\\d ( £w4 L ^ {B{p o. R)) 



In particular, 



- « + exp (^) - «_ ^ exp (-^ 



cosh 



K h%-R 



which provides a bound of c^" -1 " \4 (p°) in terms of lower derivatives given by 

di n 3 +1) K(p°) < c\R\\dg>W4L- lB < J fl,R)) + ~H?V4L~ mP0 , R)) 



(112) 



that 



Using (112 1 for the first term on the right hand side of (111), we see for all R sufficiently small 

C l+1) ^*(P°) 



< C 



R\\di?V4L~(B( p im + ^llC l)W/ *IU-(B( P0 ,fl)) 



This inequality together with (112) shows that 

( ° \\\d^W4 L ^ iB{p o iR)) + \\di?V4 L ~ {B{po , m ] . (113) 



< 



R 



To conclude, we will show by induction that (113) implies 

C,?n n 



0<?>W.(p°) + di n M(p°) 



< 



R' 



\\W*\\l°°(b( p °,r)) + \\V*\\l<*>(b( Po m)) 



(114) 
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Since this holds for any value of n, analyticity of W* and V* in the X3 direction then follows from 
Sterling's formula. 

By (1131, the above estimate holds for n = 1. Assume that it holds for some n. Using estimate 
([113| we find, 



€: +1) W*(p ) + ^ +1) K(p°) 



< 



C {n + 1) 
R 



Note that every point in B(p , is contained in a ball of radius which is contained in 



B(p°,R). Thus, applying (114) to the right hand side we obtain 

^ +1 )vM/)| + |^ +1) Kb°) 

C (n + 1) CJn 



< 



R 



( nR \r 



CZ +1 (n+l) n+1 
RP+ 1 



||W*||l«>(.B(p°,R) + ||K||l°°(B( Po ,_R) 
I|W*||l°°(B( p °,.R) + ||K||l°°(B(po,H) 



which is our induction thesis. 



A. 2 Proof of Proposition [5] 



Most properties are easily verified by inspection. Regarding cAi let us check its regularity in the 
support of VcJ" x [-L, L]. The definition of </>i being local, it suffices to look at one cell, centered in 
(0,0). We have 

d^i = d^ L c% + (1 - <) d^L (0, 0, aj 3 ) , 



and for i = 1, 2, 



1 



e 2 In (e a /(2er E )) x\ + xt 



(xidi<p L (C, x 3 ) + x 2 d 2 4> L ((, X3)) 



^cci^il + ii^iu, 

since, thanks to 



e I \nie' y - l /r e ) < C. 



The conclusion for cj>i follows. 
Regarding 02, note that 



<h(x)= X] cl)^ m ,n,e{x){(j)L{me,n£,xz) - 4>l{x)) 

(m,n)£l c 

Again, the definition of 02 is local. In the cell [—e/2, e/2] 3 , we have 



M <CV||V0 L |[ 



L°°(f2) 



for r < er F 



^2 = for e a < 2r, 



with r = y/ x\ + x\. Thanks to ([3]), this yields the global bound 

\U <Ce° +2 \\VM\ L ~(n)- 
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A. 3 Proof of Lemma [9] 

Without loss of generality, we can assume that the centre of D m ^ n ^ e is the origin. 
Note that from ( 88 1 we know that 

||We|U°°(a,\£> e ) < SUp|u m ,„| =: M . 

Remember that two linearly independent radial solutions of 

-Au + Xu = 

for \x\ > are I (v%) and K (VX-^j, the modified Bessel functions of the first and second kind. 
They are defined by 

00 

r cosh t 



K (r) = / e~ rcostlt dt , (115) 



OO/9 



n—0 V / 



see, e.g. [32] ■ Introduce 



where a and /3 are chosen so that 

ip(er e ) = u m . n , ip(e/2) = M. 

Thanks to the maximum principle, u s {x) < ip(x). Since the two functions agree on dD mtn £l the 
normal derivative of tp at |x| = sr e gives an upper bound for the normal derivative of u e at dD m>n y. 

nil 

sup ^<|V/( E r e )|. (117) 
dD m ,„. E on 

We thus proceed to estimate \ijj'(er e )\. The constants a and /? are given by 

u m , n K Q (yXs/2\ - MK (V\r e e\ 



m,nio (VAe/2) + M7 (VAr e e) 



P = 



It 



The proof relies on precise estimates on a and /3 for various regime of A. 

It is convenient to introduce the notation z = y\e. Let us hrst consider the case when 



z < r e 3 , 



We shall use an ad-hoc bound, easily verifiable using the Frobenius decomposition of Kq, see e.g. 
[3"2"] . For all x > 0, we have 

K (x) = I (x) (- In (|) - 7 + + e)^-^)^ , (118) 
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where t(x) — 4+x s/2^/2^ eX anc ^ satisfies | < Rq(x) < | for all x > 0. We now compute 



I (r £ z)K (lz)-K {r e z)h(\z) 
lo (\z) h (r e z) 

1 N *(§*) D A 



= In (2r e ) + In {-z + ej ^^^0 - In (r e * + e) *o (r.*) . 



t (r e z) 



Note that 



< |max(l, ln(a;))| 

< i| ln(r £ )| when x < 



Therefore we have 



with 



l,(r :)/v (1 ( l -z\ -K Q {r e z)lJ l -z\ = I (±z ) U>-Om(r ) (H E_- ) . (UO) 



\E £ \ < - for all z < r £ ^ 3 . 



We can now estimate a and (3 as follows 



la < 



+ 



K (z/2) 



I Q {r e z)K (\z)-K {r e z)h(\z) 
K (r £ z) 



M 



< M 



+ M 



I {r £ z)K (\z)-K {r £ z)I (\z) 
-ln(f)- 7 + ln(f +e)40 T) R o (z/2) 



I (r £ z)\ln(r £ )\(l + E £ ) 
-In (r £ z) - + ln(r £ z + e) T ^0^R o (r e z) 



M 

< 3— r +9 



h(r e z) h{\z) 



I (r £ z)\ln(r £ )\(l + E £ ) 

M _-,,o 
when z <r r 1 . 



Therefore 



Similarly, we obtain 



< 



Id < C 



M 



h (r £ z) 



M\I (z/2)\+M\I (r £ z)\ 



< CM. 



< C 



M 



M 



\h(r £ z)K (\z)-K {r £ z)I (\z)\- h (r £ z) | In(r e )| " \\n(r £ )\- 
We are now in position to bound V', namely 
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At t — er £ , since \/Xer £ < r £ < 1, we deduce that 1\ y\/XtJ < r s . Using the fact that for all 
x > 0, xK\(x) < 1, we obtain 

< V\CM i ; ' ; 

We have thus shown that 



< CM- 



VXer e \\n(r £ )\J " er e |In(r B )| 



W{er e )\ < C6(\,e)M for z < r £ s . 

— 1/3 

Let us now turn to the case when z > r £ .In this case, we easily show that 

1 . „ 1 



a s» M- 



leading to the bounds 



For r e < z < r_ 1 , we have 



, 1 , and p ~ u m n , , 
Mr*) K {r e z) 

\ner,)\<^CM(^ + ^\ 
\lo(kz) K (r £ z) 



Ml*) Io(|VAe) 



as e -> 0. 



On the other hand, it is easy to check that for < x < 1, 

Ki (x) < 2 1 
i^o(x) " i In (| 

We thus obtain for z = r~ Q , | < a < 1, that 



|^(er e )| <CM 1 + - -!— 

For the remaining case z > r^ 1 , note that 



< CM 



r e (l — a + e 2 ) 



Zi(r e x) < I (r £ x) < I y-x ) for x > 0, 
K x {x) < 2K Q (x) for x > 1. 



Inserting these inequalities in (121), we obtain 



From ( |117[ ), ( [120] ), ([122] and ( |123| , we conclude that 



|i/;'(er £ )| < CA/v 7 ! for z > r" 1 . 
conclude that 

< |^(er e )| <C0(A,e)M. 



sup — 
9D m „ e on 



(120) 



(121) 



= C6>(A,e)Af. (122) 



(123) 



The lower bound for §^ is obtained by similar arguments. Since u £ is constant on dD m n ^ £ , this 
concludes the proof. 



52 



A. 4 Proof of Lemma [lOl 

We continue to use M = sup m „ |u,„ irl |. Fix a multi-index J g {(0, 1), (1, 0)} and split \7 J u £ (x) — 
I({x) + l£(x) + l£(x) as in @. 
We first show that 



[Vtt e ]c(u.«) < CM for r > re^ 1 "^/ 2 . 



(124) 



Recall that we have shown using (96) and (97) with (3 = that 

|V/i 7 (rr)| < CM and |V/ 2 7 (x)| < 



Fix x and z in uj^. and € (0, 1). Then the above estimates implies that 
\l{(x)-l{(z)\ + \l£(x)-l£(z)\ 



< CM 



To establish (124), it remains to bound 

\l£(x)-l£(z)\ 



A(x, z) 



\x — z\ l 



= £ 

(m,n)e/ e 



d ^J y ) [V^(x - y) - V^(z - y)] 



dD„ 



dn 



For A > — , thanks to (|T00b, we have 



| a; — z| l 



1/2 



d<r(y). 



< c\vi£\ = c\ii +1 \ < cm^- 5 [^ + 1). 

which is sufficient provided r > ke 2 for example. We henceforth assume that A < We hrst bound 

the derivatives of//. For £ G uj e t1 we use Lemma[9]and count contributions from "rings" of inclusions, 
distinguishing near and far contributions as before, to get 



Ce- 1 

|V/ 3 J (e)| < C J2 e(X,e)M(j + l)\\K ( 2 \VX(j+T)e 
<CM8(\,e)er s \ i^ 2) (\/AT£)| 

Ce- 1 

+ CM9{X,e)er s ^ MjV\ K^ 2) (VXje) 

i=i 

Cs- 



<CM6(\,e)sr E — ^ + — 
I re 2 7 e 

3=1 

1 , lloge 



< CM6{\ 1 e)er e -— + 



Noting that 0(A, e) < ^- by (91 ) and our restriction on A, we thus have 



|W 3 J (£)| < CM K + Uoge 



(125) 
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From ( 125 ) we deduce in particular that 

A(x,z) <CMe 1 -" +|loge|) for \x - z\ < lOe. 



(126) 



We turn to bounding A when \x — z\ > lOe. Pick x and z in uj% such that \x — x\ < 2e, \z— z\ < 2e, 
dist (x, D e ) > g/10, dist (z, D e ) > e/10 and \x - z\ < \x - z\. Then by fll26]), 



A(ar,z) < A(ar, 5) + z ) + A(z,z) < + CM e 1 ^ [\ + \ loge|V 

Thus, provided r > keVj can focus on case when dist (x,D e ) > e/10 and dist (z, D e ) > e/10. 
Split A(x, z) — Ai(x, z) + A 2 {x, z) where 



A 1 {x,z)= 22 



A 2 (x,z)= 22 



(m,n)£:I £ (z,x) 



gfi, [VMx-y)-V^(z-y)] 

0tt e [V^(x - y) - Vj$(z - y)] 
dD mne on \x-z\ v 



and 



I e (x, z) = {(m, n) € J E : dist (a;, f m , n , £ ) < dist (z, D TO , n , e )} . 

In the sequel we bound A\(x, z) independently of x and z. Switching the role of x and z, we can 
therefore use the same bound for A 2 (x, z). 

Using Lemma [9] and the expression for 0(X,e) for A < i, we have 



\Ai(x,z)\< CMe 2 



sup 



[Vi$(x -y)- Vt$(z - y)] 



(127) 



It is convenient to introduce 



s.s'X | S — S | 



We will use the following inequality, which is easily proved using the monotonicity properties of Kq 
and Kq 2 \ the derivative of 



K { ^\t)<-^iort>Q 



C 



(128) 



To bound A\(x, z), we proceed as before by counting contribution from inclusions located in the 
rings centred at x with radii (j + l/10)e and (j + 1 + l/10)e. For example, consider an inclusion 
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D m ,n,e which is closer to x than to z and lying in the above j-th ring. For y € dD m rit£ , we estimate 
\V J x $(x-y)-V J x <P(z-y)\ 



1 



| a; - z\ v 



VXKM(SX\ X - y\) _ VAK£\JX\ Z - y\) ^-f 



< 



(x-y)-J (z-y)-J 



\x-y\ \z~y\ 
(!)/ 



V\K^(VX \x - y\) - VXK^iVX \z - y\) 



<^X\K^(V\(j + l/10)e)\ 



C 



((j + l/10)e)» 



1^ — 2/1 

+ ^ 1+iy ^ 1+ ^(VA(j + l/10) £ ) 



< 



C 



((i + 1/10) 



Inserting this estimate in (127) and summing over all inclusions we end up with 
A l {x,z)<Yj-^—i— T —<CMe 1 - v Y —<CM. 

Likewise 



A 2 (x,z) <CM. 



We have thus established (124). 



.0 



To establish the result, we note that, by Lemma \7\ we have, for < (3 < \ and r > Ke 2 i'+ 2 « , 

C 



This implies that 



|VG e |U~ K) < — ~ A/. 



|Vm e (x)-Vu £ (z)| < C M 



A/3 |:r-z| l 



On the other hand, for some small S > to be determined, (124) implies that for r > ns^ that 

\Vu £ (x) ~VU e (z)\ „a-6-u < r m u .l 1 -- 5 -^ 
i ~ S [Vu e Jc*i-,S(- wE ) \x — z\ < O M |x — z\ 



\x - z\" 

It follows that, provided < 1 — S — v, 



\Vu e {x)-Vu e {z)\ < f C M 
\P \x- 
C 



M provided r > n e V 2 1+23 J . 



g(i_a_„) 

A l ~ s 



Picking /3 = ^Ht- and S = \— ff = il 9l ^„ , we get the assertion, provided ^ < 2/3. 

1 l_|_2p i-|-zp /y 
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